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Abstract. Whether integrable, partially integrable or nonintegrable, nonlinear 
ir\ '' partial differential equations (PDEs) can be handled from scratch with essentially 

r^ , the same toolbox, when one looks for analytic solutions in closed form. The basic 

' tool is the appropriate use of the singularities of the solutions, and this can be done 

2 j ' without knowing these solutions in advance. Since the elaboration of the singular 

manifold method by Weiss ei ai, many improvements have been made. After some 
basic recalls, we give an interpretation of the method allowing us to understand 
why and how it works. Next, we present the state of the art of this powerful tech- 
nique, trying as much as possible to make it a (computerizable) algorithm. Finally, 
PlJ ' we apply it to various PDEs in 1 + 1 dimensions, mostly taken from physics, some 

of them chaotic : sine-Gordon, Boussinesq, Sawada-Kotera, Kaup-Kupershmidt, 
complex Ginzburg-Landau, Kuramoto-Sivashinsky, etc. 
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2 R. Conte 

1 Introduction 

Our interest is to find explicitly the "macroscopic" quantities which ma- 
terialize the integrability of a given nonlinear differential equation, such as 
particular solutions or first integrals. We mainly handle partial differential 
equations (PDEs), although some examples are taken from ordinary differen- 
tial equations (ODEs). Indeed, the methods described in these lectures apply 
equally to both cases. 

These methods are based on the a priori study of the singularities of the 
solutions. The reader is assumed to possess a basic knowledge of the singu- 
larities of nonlinear ordinary differential equations, the Painleve property for 
ODEs and the Painleve test. All this prerequisite material is well presented in 
a book by Hille |Q while Cargese lecture notes [g6| contain a detailed exposi- 
tion of the methods, including the Painleve test for ODEs. Many applications 



are given in a review |10S] 



As a general bibliography on the subject of these lectures, we recommend 
Cargese lecture notes [M and a shorter subset of these with emphasis on the 
various so-called truncations p4[ . 

Throughout the text, we exclude linear equations, unless explicitly stated. 

2 Various levels of integrability for PDEs, definitions 

In this section, we review the required definitions (exact solution, Backlund 
transformation. Lax pair, Darboux transformation, etc). 

The most important point is the global nature of the information which is 
looked for. The existence theorem of Cauchy (for ODEs) or Cauchy-Kowalevski 
(for PDEs) is of no help for this purpose. Indeed, it only states a local prop- 
erty and says nothing on what happens outside the disk of definition of the 
Taylor series. Therefore it cannot distinguish between chaotic equations and 
integrable ones. 

Still from this point of view, Laurent series are not better than Taylor 
series. For instance, the Bianchi IX cosmological model is a six-dimensional 
dynamical system 

cr^ {hog A)" = A^ -{B-Cf, and cyclically, cr^ = 1, (1) 



which is undoubtedly chaotic |107|. Despite the existence of the Laurent series 

n 

A/o- = X'^ + a2X + O(x^), X^T-Ti, 

B/n = b^X + biX^+0{x^), (2) 

which depends on six independent arbitrary coefficients, (ti, 5o, cq, 6i, ci, 02), 
a wrong statement would be to conclude to the absence of chaos. 
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This leads us to the definition of the first one of several needed global 
mathematical objects. 

Definition 1. One calls exact solution of a nonlinear PDE any solution 
defined in the whole domain of definition of the PDE and which is given in 
closed form, i.e. as a finite expression. 

The opposite of an exact solution is of course not a wrong solution, but 
what Painleve calls "une solution illusoire" , such as the above mentioned 
series. 

Note that a multivalued expression is not excluded. From this definition, 
an exact solution is global, as opposed to local. This generically excludes 
series or infinite products, unless their domain of validity can be made the 
full domain of definition, like for linear ODEs. 

Example 2. The Kuramoto-Sivashinsky (KS) equation 

Ut + UU.^ + ^lU.^x + VU^xxx =0, i^ 7^ 0, (3) 

describes, for instance, the fiuctuation of the position of a fiame front, or the 
motion of a fluid going down a vertical wall, or a spatially uniform oscillating 
chemical reaction in a homogeneous medium (see Ref. |7q] for a review) , and 
it is well known for its chaotic behaviour. An exact solution is the solitary 
wave of Kuramoto and Tsuzuki p39| in which the wavevector k is fixed 

/ k fc\'^/60 \ k k 

u = UOiy - tanh -M + — ^ - SOi'k'^ - tanh -^ + c, 
\ 2 ^ / \ 1"^ / 

^ = x-ct-xo, k^ = — — or - — -, (4) 

which depends on two arbitrary constants (c, xq). On the contrary, the Lau- 
rent series 

u = 120.r^ + ^Mr ^ + c - ^M^C + "6^^ + 0{i% (5) 

which depends on three arbitrary constants {c,Xq,Uq), is not an exact solu- 
tion, since no closed form expression is yet known for the sum of this series. 

There exists a powerful tool to build exact solutions, this is the Backlund 
transformation. For simplicity, but this is not a restriction, we give the basic 
definitions for a PDE defined as a single scalar equation for one dependent 
variable u and two independent variables {x,t). 



Definition 3. (Refs. |41| vol. Ill chap. XII, |M) A Backlund transfor- 



mation (BT) between two given PDEs 

El {u, x,t) = Q, E2 [U, A, T) = (6) 
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is a pair of relations 



Fj{u,x,t,U,X,T) = 0, j = 1,2 (7) 



with some transformation between {x,t) and (X,T), in which Fj depends on 
the derivatives oiu{x,t) and U(X,T), such that the ehmination of u (resp. U) 
between (-^1,^2) impUes E2{U,X,T) = (resp. Ei{u,x,t) = 0). The BT is 
called the auto-BT or the hetero-BT according as the two PDEs are the 
same or not. 

Example 4- The sine-Gordon equation (we identify sine-Gordon and sinh- 
Gordon since an afhne transformation on u does not change the integrability 
nor the singularity structure) 

sine-Gordon : _E(w) = w^t + 2asinhu = (8) 

admits the auto-BT 

u-U 
(u + [/)a; + 4Asinh— — =0, (9) 

(«-C/).-f sinh^^O, (10) 

in which A is an arbitrary complex constant, called the Bdcklund parameter. 

Given the obvious solution U — Q (called vacuum), the two equations 
(Pl)~(M) ^I'S Riccati ODEs with constant coefficients for the unknown e"'^, 

{e^'% = A(l - (e"/2)2)^ (11) 

(e"/2), = -a(l-(e"/2)2)/(2A), (12) 

therefore their general solution is known in closed form 

e"/2 = tanhe, 9 = (\x ~ -^t - zq) , (13) 

with (A, Zq) arbitrary. This solution is called the one-soliton solution, it is 
also written as 

tanh(w/4) = -e~^^, Ux = 4Asech26', ut = -2aA"^ sech26'. (14) 

By iteration, this procedure gives rise to the A^-soliton solution 170,01, an 
exact solution depending on 2N arbitrary complex constants {N values of 
the Backlund parameter A, N values of the shift zq), with N an arbitrary 
positive integer. A remarkable feature of the SG-equation, due to the fact 
that at least one of the two ODEs (|9|)-([l0|) is of order one, is that this N- 
soliton can be obtained from N different copies of the one-soliton by a simple 
algebraic operation, i.e. without integration (see Musette's lecture [p5[). 
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Example 5. The Liouville equation 

Liouville: E{u) = Uxt + ae" = (15) 

admits two BTs. The first one 

(w-v), = aAe("+")/^ (16) 

(w + z;)t = -2A-ie("-'')/2^ (17) 

is a BT to a linearizable equation called the d'Alembert equation 

d'Alembert: E{v) = v^t = 0. (18) 

The second one is an auto-BT 

u -U 
(u + C/), = -4Asinh^— , (19) 

(u-C/)t = A-iae("+^)/2. (20) 

The first of these two BTs allows one to obtain the general solution of the 
nonlinear Liouville equation, see Section |7|. 

This ideal situation (generation of the general solution) is exceptional 
and the generic case is the generation of particular solutions only, as in the 
sine- Gordon example. 

The importance of the BT is such that it is often taken as a definition of 
integrability. 

Definition 6. A PDE in N independent variables is integrable if at least 
one of the following properties holds. 

1. It is linearizable. 

2. For A^ > 1, it possesses an auto-BT which, if A^ = 2, depends on an 
arbitrary complex constant, the Backlund parameter. 

3. It possesses a hetero-BT to another integrable PDE. 

Although partially integrable and nonintegrable equations, i.e. the ma- 
jority of physical equations, admit no BT, they retain part of the properties 
of (fully) integrable PDEs, and this is why the methods presented in these 
lectures apply to both cases as well. For instance, the KS equation admits the 
vacuum solution u — and, in Section 0, an iteration will be built leading 
from u = to the solitary wave (0); the nonintegrability manifests itself in 
the finite number of times this iteration provides a new result (A^ = 1 for the 
KS equation, and one cannot go beyond (0) pOl). 



For various applications of the BT, see Ref. [104|. 

When a PDE has some good reasons to possess such features, such as the 
reasons developed in Section 0, we want to find the BT if it exists, since this 
is a generator of exact solutions, or a degenerate form of the BT if the BT 
does not exist, and we want to do it by singularity analysis only. 

Before proceeding, we need to define some other elements of integrability. 
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Definition 7. Given a PDE, a Lax pair is a system of two linear differential 
operators 

Laxpair: Li(t/,A), L2((7,A), (21) 

depending on a solution U of the PDE and, in the 1 + l-dimcnsional case, on 
an arbitrary constant A, called the spectral parameter, with the property 
that the vanishing of the commutator [Li,Z/2] is equivalent to the vanishing 
of the PDE E{U) = 0. 

A Lax pair can be represented in several, equivalent ways. 
The Lax representation |72|| is a pair of linear operators (L, P) (scalar or 
matrix) defined by 

Li=L-X, L2 = dt-P, Lii, = 0, 121^ = 0, Xt = 0, (22) 

in which the elimination of the scalar A yields 

Lt^[P,L], (23) 

i.e. , thanks to the isospectral condition At = 0, a time evolution analogous 
to the one in Hamiltonian dynamics. 

The zero-curvature representation is a pair (i, M) of linear operators 
independent of {dx,dt) 

Li^dx- L, L2 = dt- M, Liip = 0, LaV' = 0, 

[dx ~ L,dt- M] ^ Lt- Mcc + LM - ML = 0. (24) 

The common order N of the matrices is called the order of the Lax pair. 

The projective Riccati representation is a first order system of 27V — 2 
Riccati equations in the unknowns V'j/V'ijJ — 2,...,A'^, equivalent to the 
zero-curvature representation (|24|). 

The scalar representation is a pair of scalar linear PDEs, one of them of 
order higher than one, 

Liij = 0, L21P = 0, 

X=[Li,L2] = 0. (25) 

In 1 + 1-dimensions, one of the PDEs can be made an ODE (i.e. involving 
only X- or t-derivatives), in which case the order of this ODE is called the 
order of the Lax pair. 

The string representation or Sato representation ||64| 

[P,Q]^1. (26) 

This very elegant representation, reminiscent of Hamiltonian dynamics, uses 
the Sato definition of a microdijferential operator (a differential operator 
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with positive and negative powers of the differential operator d) and of its 
differential part denoted ()+ (the subset of its nonnegative powers), e.g. 

Q = dl-u, (27) 

L^Q^'^, (28) 

{L^)^^dl--{?,/A){u,d^}^ (29) 

{L")^ ^dl~ (5/4){u, dl} + (5/16){3u2 + u^^, 9J, (30) 

in which {a, 6} denotes the anticommutator ab+ba. See Ref. [E2| for a tutorial 
presentation. 

Example 8. The sine-Gordon equation (^ admits the zero-curvature repre- 
sentation 



(9,-M)(^Ji)=0, M=-(a/2)A-i(^ 



e^ 
e-^ 



(32) 



equivalent to the Riccati representation, with y — 'ipi/'4'2, 



2 



Vx ^ >^ + Uxy- Xy, (33) 

y* = -^A-ie^ + ^A-ie-V. (34) 

Example 9. The matrix nonlinear Schrodinger equation 

iQt + {b/a)Q^^ - 2abQRQ = 0, -iRt + {b/a)R^^ - 2abRQR = 0, (35) 

in which (Q, i?) are rectangular matrices of respective orders (m, n) and 
(n,m), and {i,a,b) constants, admits the zero-curvature representation ([ [77| 
Eq. (5)) 

{dx - L)^ = 0, {dt - M)V = 0, (36) 

L = aP + AG, Af = i-aGP^ + GP^ + 2AP + {2/a)X'^G)b/i, (37) 

in which A is the spectral parameter, P and G matrices of order m + n defined 
as 

The matrix G characterizes the internal symmetry group GL(to, C)(8'GL(n, C). 
The lowest values 

m = 1,71= 1, Q = (u), i?= ([/), (39) 



define the AKNS system (Section 9J_), whose reduction U = u is the usual 
scalar nonlinear Schrodinger equation. 
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Example 10. The 2 + 1-diniensional Ito equation 163] 

E[u) = {u^xxt + Qa'^UxtUxx + aiUtt + a2U^t + a^iU^^ + a^Uty) ^ = (40) 

has a Lax pair whose scalar representation is 

Li = dl + aidt + (aa + 6a-i[/^^)<9^ + a^dy - X (41) 

L2 = d^dt - ^x^^ + (y + 2«"'C/,t) (42) 

a[Li,L2]-2^([/) + 6C/,,,L2- (43) 

In the 2 + 1-dimensional case 04 7^ 0, the parameter A can be set to by the 
change ip 1-^ 'tpe^^ . This is the reason of the precision at the end of item 2 in 
definition H. This pair has the order four in the generic case oi 7^ 0, although 
neither Li nor L2 has such an order. 

Example 11. The string representation of the Lax pair of the derivative of 
the first Painleve equation is 

[P, Q\ = [{{dl ~ uf)^ .dl~u]^ -(l/4)«... + {'i/A)uu.,, = 1. (44) 

Exam,ple 12. The Sato representation of the Lax pair for the whole Korteweg- 
de Vries hierarchy is 

9t„L=[(L2™-i)^,i], L = gi/2^ Q^dl~u,m^\,?,,^,... (45) 

From the singularity point of view, the Riccati representation is the most 
suitable, as will be seen. 

The last main definition we need is the Darboux transformation (DT). 
The working definition given below is very simplified (this is an involution) 
as compared to the one of Darboux [|40| , but it is sufficient for our purpose. 
The full definition is given in Musette's lecture |85|| . 

Definition 13. Given a PDE, a Darboux transformation is a transfor- 
mation between two solutions (u, U) of the PDE 

DT : u^^Vf Logr/ + U (46) 

/ 
linking their difference to a finite number of linear differential operators Vf 
(/ like family) acting on the logarithm of functions ry. 

In the definition (|46|), it is important to note that, despite the notation, 
each function t/ is in fact the ratio of the "tau-function" of u by that of U . 

Lax pairs, Backlund and Darboux transformations are not independent. 
In order to exhibit their interrelation, one needs an additional information, 
namely the link 

V/ : VfhogTf^Ffii;), (47) 

which most often is the identity t = tp, between the functions ry and the 
function Tp in the definition of a scalar Lax pair. 
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Example 14- The (integrable) sine-Gordon equation (H) admits the Darboux 
transformation 

w = t/-2(Logri-LogT2), (48) 

in which (ti, T2) is a sohition {ipi, ^2) of the system (pi|)-(p2|). 

Then its BT (||)-(|l^) is the result of the ehmination fl^ o f ti/t2 between 
the DT (14^ ) and the Riccati form of the Lax pair (|3i^)-(|34|), with the cor- 
respondence Tf = tpf,f = 1,2. This elimination reduces to the substitution 
y — e~^"^^'>^^ in the Riccati system (|3it)-(p4|), and this is one of the advan- 
tages of the Riccati representation. Therefore the Backlund parameter and 
the spectral parameter are identical notions. 

Example 15. The (nonintegrable) Kuramoto-Sivashinsky equation admits the 
degenerate Darboux transformation 

u^U + {mvdl + {m/l9)^jLd^) Log T, (49) 

in which U = c (vacuum) and r is the general solution V' of the linear system 
(a degenerate second order scalar Lax pair) 

LiiP = [dl - /cV4)V' = 0, (50) 

L2^ = (dt + cd^)^ = 0, (51) 

[ii,i2] = 0. (52) 

The solution u defined by (09) is then the solitary wave (0), and this is 
a much simpler way to write it, because the logarithmic derivatives in (Eof) 
take account of the whole nonlinearity. 

Since, roughly speaking, the BT is equivalent to the couple (DT, Lax 
pair), one can rephrase as follows the iteration to generate new solutions. Let 
us symbolically denote 

E{u) = the PDE, 

Lax(^, A, C/) = a scalar Lax pair, 

F the link <^ VLogT = F{ip) from ■0 to r, 

u = Darboux([/, r) the Darboux transformation. 

The iteration is the following, see e.g. [p3. 

1. (initialization) Choose uo = a particular solution of E{u) = 0; set ji = 1; 
perform the following loop until some maximal value of n; 

2. (start of loop) Choose A„ — a particular complex constant; 

3. Compute, by integration, a particular solution ijjn of the linear system 
Lax(V', A„,M„_i) = 0; 

4. Compute, without integration, I?LogT„ = F{ipn)] 

5. Compute, without integration, u„ = Darboux(M„_i, t„); 

6. (end of loop) Set n — n + 1. 

Depending on the choice of A„ at step 2, and of ^n at step 3, one can 
generate either the iV-soliton solution, or solutions rational in {x,t), or a 
mixture of such solutions. 
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3 Importance of the singularities : a brief survey of the 
theory of Painleve 

A classical theorem states that a function of one complex variable without 
any singularity in the analytic plane (i.e. the complex plane compactifled by 
addition of the unique point at infinity) is a constant. Therefore a function 
with singularities is characterized, as shown by Mittag-Leffler, by the knowl- 
edge of its singularities in the analytic plane. Similarly, if u satisfies an ODE 
or a PDE, the structure of singularities of the general solution characterizes 
the level of integrability of the equation. This is the basis of the theory of 
the (explicit) integration of nonlinear ODEs built by Painleve, which we only 
briefly introduce here [for a detailed introduction, see Cargese lecture notes : 
Ref. (2| for ODEs, Ref. [|4| for PDEs]. 

To integrate an ODE is to acquire a global knowledge of its general so- 
lution, not only the local knowledge ensured by the existence theorem of 
Cauchy. So, the most demanding possible definition for the "integrability" of 
an ODE is the single valuedness of its general solution, so as to adapt this 
solution to any kind of initial conditions. Since even linear equations may 
fail to have this property, e.g. 2xu' -I- m = 0, u = cx~^''^, a more reasonable 
definition is the following one. 

Definition 16. The Painleve property (PP) of an ODE is the uniformizabil- 
ity of its general solution. 

In the above example, the uniformization is achieved by the change of in- 
dependent variable x = X'^. This definition is equivalent to the more familiar 
one. 

Definition 17. The Painleve property (PP) of an ODE is the absence of 
movable critical singularities in its general solution. 

Definition 18. The Painleve property (PP) of a PDE is its integrability 
(Definition 0) and the absence of movable critical singularities in its general 
solution. 

Let us recall that a singularity is said movable (as opposed to fixed) if 
its location depends on the initial conditions, and critical if multivaluedness 
takes place around it. Indeed, out of the four configurations of singularities 
(critical or noncritical) and (fixed or movable), only the configuration (critical 
and movable) prevents uniformizability : one does not know where to put the 
cut since the point is movable. 

Wrong definitions of the PP, alas repeatedly published, consist in replacing 
in the definition "movable critical singularities" by "movable singularities 
other than poles" , or "its general solution" by "all its solutions" . Even worse 
definitions only refer to Laurent series. See Ref. ^^, Section 2.6, for the 
arguments of Painleve himself. 
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The mathematicians hke Painleve want to integrate whole classes of ODEs 
(e.g. second order algebraic ODEs). We will only use their methods for a given 
ODE or PDE, with the aim of deriving the elements of integrability described 
in Section (exact solutions, . . . ). This Painleve analysis is twofold ("double 
methode", says Painleve). 

1 . Build necessary conditions for an ODE or a PDE to have the PP (this is 
called the Painleve test). 

2. When all these conditions are satisfied, or at least some of them, find 
the global elements of integrability. In the integrable case this is achieved 
either (ODE case) by explicitly integrating or (PDE case) by finding 
an auto-BT (like equations (^)-(foh for sine-Gordon) or a BT towards 



another PDE with the PP (like (h6h-(17) between the d'Alembert and 
Liouville equations). In the partially integrable case, only degenerate 
forms of the above can be expected, as described in Section 0. 

4 The Painleve test for PDEs in its invariant version 

When the PDE reduces to an ODE, the Painleve test (for shortness we will 
simply say the test) reduces by construction to the test for ODEs, presented 
in detail elsewhere [^ and assumed known here. 

We will skip those steps of the test which are the same for ODEs and 
for PDEs (e.g., diophantine conditions that all the leading powers and all 
the Fuchs indices be integral) , and we will concentrate on the features which 
are specific to PDEs, namely the description of the movable singularities, the 
optimal choice of the expansion variable for the Laurent series, the advantage 
of the homographic invariance. 

4.1 Singular manifold variable cp and expansion variable % 

Consider a nonlinear PDE 

E{u,x,t,...)^0. (53) 

To test movable singularities for multivaluedness without integrating, 
which is the essence of the test, one must first describe them, then, among 
other steps, check the existence near each movable singularity of a Laurent 
series which represents the general solution. 

For PDEs, the singularities are not isolated in the space of the independent 
variables {x, t, . . .), but they lay on a codimension one manifold 

ipix,t,...)-ipo^O, (54) 

in which the singular manifold variable (p is an arbitrary function of the 
independent variables and ipo an arbitrary movable constant. Even in the 
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ODE case, the movable singularity can be defined as >p{x) — ipo — 0, since the 
implicit functions theorem allows this to be locally inverted to x — xo =0; the 
arbitrary function ip thus introduced may then be used to construct exact 
solutions which would be impossible to find with the restriction (p{x) — x 

One must then define from (p — ifQ an expansion variable x for the Laurent 
series, for there is no reason to confuse the roles of the singular manifold 
variable and the expansion variable. Two requirements must be respected: 
firstly, X must vanish as iy9 — (/jq when ip —> ipo; secondly, the structure of 
singularities in the if complex plane must be in a one-to-one correspondence 
with that in the x complex plane, so x must be a homographic transform of 
if — ipo (with coefficients depending on the derivatives of ip) . 

The Laurent series for u and E involved in the Kowalevski-Gambier part 
of the test are defined as 



-t-CXD 



+ OC 



" = E "jX^'^^ -peM, E^Y. ^^■^'^'' -« ^ ^* (55) 

with coefficients Uj , Ej independent of x s-ud only depending on the deriva- 
tives of ip. 

To illustrate our point, let us take as an example the Korteweg-de Vries 
equation 



E = but 



{6/a)uUx — 



(56) 



(this is one of the very rare locations where this equation can be taken as an 
example; indeed, usually, things work so nicely for KdV that it is hazardous 

to draw general conclusions from its single study). 

The choice x = V ~ 'Po originally made by Weiss et al. [ 119 makes the 
coefficients Uj , Ej invariant under the two-parameter group of translations 
ip i-^ ip + b' , with b' an arbitrary complex constant and therefore they only 
depend on the differential invariant grad ip of this group and its derivatives : 



2a<^'x"'- 



-1 '-Pt 

-2aipxxX +ab—— 
bipx 



3 ipx 2 



Px 



fx 



+ 0(x), X^V-Vo- 

. (^^) 
There exists a choice of x for which the coefficients exhibit the highest 

invariance and therefore are the shortest possible (all details are in Section 

6.4 of Ref. [pel), this best choice is fR3 



X = 



ip-ipo 



VXX , N 

2(y5a; 



•Px 



ip-ipQ 



2(/?a; 



, "fx^^, 



(58) 



in which x denotes one of the independent variables whose component of 
grad ip does not vanish. The expansion coefhcients Uj , Ej are then invariant 



Exact solutions by singularity analysis 19 February 2000 13 

under the six-parameter group of honiographic transformations 

^^T^' a'd'-b'c'^O, (59) 

c(p + a' 

in which a',b',c',d' are arbitrary complex constants. Accordingly, these co- 
efficients only depend on the following elementary differential invariants and 
their derivatives: the Schwarzian 

S = {ip; x} = — oi-—] , (60) 

iPx ^ \ fx J 

and one other invariant per independent variable t,y, . . . 

C = -ift/ifx, K = -ipy/(px, ... (61) 

The reason for the minus sign in the definition of C is that, under the travel- 
ling wave reduction £^ — x~ ct, the variable C becomes the constant c. These 
two invariants are linked by the cross-derivative condition 

X = {i(pxxx)t - {vt)xxx)/Vx = St+ Cxxx + "^.CxS + CSx = 0, (62) 

identically satisfied in terms of (p. 

For our KdV example, the final Laurent series, as compared with the 
initial one (p7|), is remarkably simple : 



u = 2ax-' -ab- + ^-~ 2a{hC - S)xX + O(x'), X = (58). (63) 

6 3 

For the practical computation of (uj,Ej) as functions of {S,C) only, 
i.e. what is called the invariant Painleve analysis, the above explicit ex- 
pressions of (S, C, x) in terms of ip are not required, the variable ip com- 
pletely disappears, and the only necessary information is the gradient of the 
expansion variable x defined by Eq. (p8[). This gradient is a polynomial of 
degree two in x (this is a property of homographic transformations), whose 
coefficients only depend on S, C: 

Xx = l + |x', (64) 

Xt^-C + CxX-\{CS + Cxx)x^. (65) 



The above choice (58) of x which generates homographically invariant 
coefficients is the simplest one, but it is only particular. The general solution 
to the above two requirements which also generates homographically invariant 
coefficients is defined by an affine transformation on the inverse of x 186|| 

Y-^ = B{x'^ + A), B^ 0. (66) 
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Since a homography conserves the Riccati nature of an ODE, the variable Y 
satisfies a Riccati system, easily deduced from the canonical one (p4)-(p5h 



satisfied by x, see (115)-(116) 



A frequent worry is : is there any restriction (or advantage, or inconve- 
nient) to perform the test with x or Y rather than with ip — ipo? The precise 
answer is : the three Laurent series are equivalent (their set of coefficients 
are in a one-to-one correspondence, only their radii of convergence arc differ- 
ent). As a consequence, the Painleve test, which involves the infinite series, 
is insensitive to the choice, and the costless choice (the one which minimizes 
the computations) is undoubtedly x defined by its gradient (|6^)~(|65|) (to per- 
form the test, one can even set, following Kruskal [p3[, S ~ 0,Cx ^ 0). If the 
same question were asked not about the test but about the second stage of 
Painleve analysis as formulated at the end of Section H, the answer would be 



quite different, and it is given in Section 6.1. 

Finally, let us mention a useful technical simplification. From its definition 
(p8|), the variable x~^ is a logarithmic derivative, so the system (|64|)-(|65|) can 
be integrated once 

'F ^ {^ - ^o)^-'/^ (67) 

(Logif), = X"\ (68) 

{Log^)t = -Cx-' + lCx. (69) 

This feature helps to process PDEs which can be defined in either conservative 
or potential form when the conservative field u has a simple pole, such as the 
Burgers equation 

E{u) = but + iu^/a + ux)x = 0, (70) 

F{v) = bvt + vl/a + «,, + G{t) =0, u = vx, E ^ F,. (71) 

Despite its (unique) logarithmic term, the ^'-series for v 

V = a Log ^ + vo + (2i;o,i — abC)x 

+ (FM - aS/2 + abC.,/2)x^ + O(x'), (72) 

in which vq is arbitrary, is "shorter" than the Laurent series for u 

u = aX"^ + {ab/2)C + U2X 

+ [{a/4){b\Ct + CCx) + 2bC ~Sx~ U2,x)] x' + O(x'), (73) 

in which U2 is arbitrary, and the resulting series for F{v), which is not a -0- 
series but a Laurent series, is much shorter than the Laurent series for E{u). 



See Section 7.3 for an application. 
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4.2 The WTC part of the Painleve test for PDEs 

As mentioned at the beginning of Section H, we do not give here all the 
detailed steps of the test nor all the necessary conditions which it generates 
(this is done in Section 6.6 of Ref. p6|). We mainly state the notation to be 
extensively used throughout next sections. 



The WTC part |119| of the full test, when rephrased in the equivalent 
invariant formalism |23| , consists in checking the existence of all Laurent series 
(p5[) ab le to represent the general solution, maybe after suitable perturbations 
[ESpS] not describe here. 

The gradient of the expansion variable x is given by (|64|)-(|65|), with the 
cross-derivative condition (p3). This condition may be used to eliminate, 
depending on the PDE, either derivatives Srax,nt, with n > 1, or derivatives 
Cmx.nt, with m > 3, and all equations later written are already simplified in 
either way. 

The first step is to find all the admissible values {p,uo) which define the 
leading term of the series for u. Such an admissible couple is called a family 
of movable singularities (the term branch should be avoided for the confusion 
which it induces with branching, i.e. multivaluedness). 

The recurrence relation for the next coefficients Uj, after replacement of 

Vj > 1 : Ej = P{j)u, + Q,{{ui \l<j})^0 (74) 

depends linearly on Uj and nonlinearly on the previously computed coeffi- 
cients Ul. 

The second step is to compute the indicial equation 

P{i) = (75) 

(a determinant in the multidimensional case of a system of PDEs) . Its roots 
are called the Fuchs indices of the family because they are indeed the char- 
acteristic indices of a linear differential equation near a Fuchsian singularity 
(the name resonances sometimes given to these indices refers to no resonance 
phenomenon and should also be avoided). One then requires that all indices 
be integral and obey a rank condition which, for a single PDE, reduces to 
the condition that all indices be distinct. The value i = — 1 is always a Fuchs 
index. 

The third and last step is to require that, for any admissible family and 
any Fuchs index i (a signed integer), the no-logarithm condition 

Vi e Z, P{i) = : Q^^Q (76) 

holds true, so that the coefficient Ui is an arbitrary function of the indepen- 
dent variables. In the multidimensional case, this is the condition of orthog- 
onality between the vector Qj and the adjoint of the linear operator P(j). 
Whenever there exist negative integers in addition to the ever present value 
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— 1 counted with multiplicity one, the condition (um can only be tested by a 
perturbation p8| . 

This ends this subset of the test which, let us insist on the terminology, 
is only aimed at building necessary conditions for the PP. 

The Laurent series for u built in this way depends on at most N arbitrary 
functions (if N denotes the differential order), namely the coefficients Ui 
introduced at the N Fuchs indices, including ip for the index —1. 

Any item uj,Ej,Qj depends, through the elementary invariants (5, C), 
on the derivatives of (p up to the order j + 1, so the dependences are as follows: 
uo = f{C), ui = f{C, C^.Ct), U2 = f{C, C,,Ct,C, C^u Ctu S),... 

Let us take an example. 

Example 19. The Kolmogorov-Petrovskii-Piskunov (KPP) equation |67| , |9^ 

E{u) = hut — Uxx + 2dr'^{u — ei){u — e2){u — 63) ~ 0, Cj distinct, (77) 
encountered in reaction-diffusion systems (an additional convection term uux 



is quite important in physical applications to prey-predator models [lOSj) 
possesses the two families {d denotes any square root of (P) 

p = -1, Uo = d, (78) 

each family has the same two indices (—1,4), and the Laurent series for each 
family reads 

u^dx^^ + {si/3-{bd/6)C) (79) 

- (d/6) ((6V6)C2 -6a2-S- bCx) x + O(x'), (80) 

with the notation 

Sl = ei -H 62 -H 63, 02 = ((62 - 63)2 + (eg - df + (d - 63)^) 7(18^2). (81) 

At index i = 4, the two no- log conditions, one for each sign of d |lq] , 

Qi = C[{bdC + si - 3ei){bdC + si - 3e2)(MC + si - Seg) 

~3b'^d^{Ct + CCx)]^0 (82) 

are not identically satisfied, so the PDE fails the test. 

It is time to define a quantity which, although useless for the test itself, 
is of first importance at the second stage of Painleve analysis, which will be 
developed in Sections o to O. This quantity is defined from the finite subset 
of nonpositive powers of the Laurent series for u. 

Definition 20. Given a family (p, uq), the singular part operator V is defined 
as 

Log ip 1-^ V Log If = ut(0) — ut(oo), (83) 
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in which the notation UT{y^o)i which emphasizes the dependence on the mov- 
able constant (po, stands for the principal part (T like truncation) of the 
Laurent series (p5|), i.e. the finite subset of its nonpositive powers 

UTi^,)=J2^jX'+''- (84) 

3=0 

For most PDEs, this operator is linear. For the Laurent series already consid- 
ered (p3), ([72|), (73), (§q), the operator is, respectively, V = — 2a9^, a, adx, ddx- 



For the Kuramoto-Sivashinsky equation (^, there exists a unique Laurent se- 
ries (55) with p — —3 (given by (||) for a particular value of x, and by the 



derivative of ( p47|) for any x)^ with a singular part operator equal to 

V = mvdl + {m/l9)^JLdx. (85) 

This is precisely the third order linear operator on the rhs of (By). 

4.3 The various ways to pass or fail the Painleve test for PDEs 

If one processes a multidimensional PDE the coefficients of which depend 
on some parameters ^, 

E(u,x;a*) = 0, (86) 

(boldface means multicomponent), the Painleve test generates the following 
output : 

1. leading order (p, uq), Fuchs indices i and singular part operator V for 
each admissible family, 

2. diophantine conditions that all singularity orders p and all Fuchs indices 
i be integral, conditions whose solution creates constraints of the type 

F(/i,C)=0, (87) 

3. no-log conditions 

Vi Vn Vuj : Qr(/i, 5, C, Uj) = 0, (88) 

arising from any integral Fuchs index i, in which n is the Fuchsian per- 
turbation order |2q| if necessary, Uj are the arbitrary coefficients Uarb 
introduced at earlier Fuchs indices j. 

In particular, the Laurent scries ( pq ) are of no use and should not be 
computed beyond the highest Fuchs integer. All this output (items 1 and 
3) is easily produced with a computer algebra program and, in all further 
examples, we will simply list these results without any more detail. 
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Strictly speaking, the answer provided by the test to the question "Has 
the PDE the PP?" is either no (at least one of the necessary conditions fails) 
or maybe (all necessary conditions are satisfied, and the PDE may possess the 
PP but this still has to be proven). It is never yes, as shown by the famous 
counterexample of Painleve (the second order ODE with the general solution 
p(ALog(cia; + 02)152,53), which therefore has the PP iff 27rjA is a period of 
the Weierstrass elliptic function p, a transcendental condition impossible to 
generate by a finite algebraic procedure). 

Now that the necessary part (i.e. the Painleve test) of Painleve analysis 
is finished, let us turn to the question of sufficiency. 

To reach our goal which is to obtain as many analytic results as possible, 
we do not adopt such a drastic point of view, but the opposite one. Instead 
of the logical and performed by the mathematician on all the necessary con- 
ditions generated by the test, we perform a logical or operation on these 
conditions. Therefore the above Painleve test must be performed to its end, 
i.e. without stopping even in case of failure of some condition, so as to collect 
all the necessary conditions. Turning to sufficiency, these conditions have to 
be examined independently in the hope of finding some global element of in- 
tegrability. An application of this point of view to the Lorenz model, a third 
order ODE, can be found in Section 6.7 of Ref. |26). 

If the PDE under study possesses a singlevalued exact solution, there must 
exist a Laurent series ( |55| ) which represents it locally. Therefore the practical 
criterium to be implemented deals with the existence of particular Laurent 
series, and the result of the test belongs to one of the following mutually 
exclusive situations. 

1. (The best situation) Success of the test, at least for some values of ^ 
selected by the test. The PDE may have the PP, and one must look for 
its BT; 

2. There exists at least one value of (/i, ip, Uarb) which ensures the existence 
of a particular Laurent series. For these values, an exact solution may 
exist; 

3. There exists at least one value of (/^,(^,Marb) enforcing some of, but not 
all, the no-log conditions of at least one particular Laurent series. Quite 
probably no exact solution exists, but there may exist a conservation law 
(a first integral for an ODE); 

4. (The worst situation) There is no value of (fi, (^,-Uarb) enforcing at least 
one of the no-log conditions of the various series. Quite probably the PDE 
is chaotic and possesses no exact solution at all. 

Examples of these various situations are, respectively : 

1. All the PDEs which have the PP (sine-Gordon, Korteweg-de Vries, . . . ), 
but also the counterexample of Painleve quoted above; 

2. The equation of Kuramoto and Sivashinsky (j^, with the particular Lau- 
rent series (ph; 
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3. The Lorenz model for b = 2a, for which the no-log condition at i = 4 is 
violated and there exists a first integral; 

4. The Rossler dynamical system for which the unique family has the never 
satisfied condition Q2 = 16 = 0. 

5 Ingredients of the "singular manifold method" 

The methods to handle the integrable and nonintegrable situations are 
the same, simply a more or less important result is obtained. 

The goal is to find a (possibly degenerate) couple (Darboux transforma- 
tion, Lax pair) in order to deduce the Backlund transformation or, if a BT 
does not exist, to generate some exact solutions. 

The full Laurent series is of no help, for this is not an exact solution 
according to the definition in Section 0. Since this is the only available piece of 
information and since a finite (closed form) expression is required to represent 
an exact solution, let us represent, and this is the idea of Weiss, Tabor and 



Carnevale [119|, an unknown exact solution u as the sum of a singular part, 
built from the finite principal part of the Laurent series (i.e. the finite number 
of terms with negative powers), and of a regular part made of one term 
denoted U. This assumption is identical to that of a Darboux transformation 
(Eq), in which nothing would be specified on U. 

This method is widely known as the singular manifold method or trunca- 
tion method because it selects the beginning of the Laurent series and discards 
("truncates") the remaining infinite part. 

Since its introduction by W TC |119| , it has been improved in many di- 
rections I 



34,101, D0|, and we present below the current status of 



the method. 



5.1 The ODE situation 

For the six ordinary differential equations (ODE) (P1)-(P6) which bear 
his name, Painleve proved the PP by showing P^,p8[ the existence of one 
(case of (PI)) or two ((P2)-(P6)) function(s) r = ti, T2, called tau-fmictions, 
linked to the general solution u by logarithmic derivatives 

(PI) : u = ViLogT (89) 

(Pn), n = 2, . . . ,6 : u^ P„(Logri - Logra) (90) 

where the operators I?„ are linear: 

Pi = -dl V2=Vi = ±d,, V3 = ±e-"a„ (91) 

V5 = ±xe-'={2a)-^^^d^, Vq = ±x{x - l)e-''{2a)-^/^d^,. (92) 

These functions ti , T2 satisfy third order nonlinear ODEs and they have the 
same kind of singularities than solutions of linear ODEs, namely they have 
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no movable singularities at all; they are entire functions for (P1)-(P5), and 
their only singularities for (P6) are the three fixed critical points (oo, 0, 1). 

ODEs cannot possess an auto-BT, since the number of independent ar- 
bitrary coefficients in a solution cannot exceed the order of the ODE. They 
can however possess a Schlesinger transformation (see definition Section O) . 



5.2 Transposition of the ODE situation to PDEs 

For PDEs, similar ideas prevail. The analogue of(p^)-(p^), with an ad- 
ditional rhs U, is now the Darboux transformation (|46|), and the scalar(s) tp 
to which the scalar(s) t are linked by (EtI) are assumed to satisfy a linear 
system, the Lax pair. 

Another interesting observation must be made. There seems to exist two 
and only two classes of integrable 1 -I- 1-dimensional PDEs, at least at the 
level of the base member of a hierarchy : those which have only one family 
of movable singularities, and those which have only pairs of families with 
opposite principal parts, similarly to the distinction between (PI) on one 
side and (P2)-(P6) on the other side. Among the 1 + 1-dimensional inte- 
grable equations, those with one family include KdV, the AKNS, Hirota- 
Satsuma and Boussinesq equations; they also include the Sawada-Kotera, 
Kaup-Kupershmidt and Tzitzcica equations because only one of their two 
families is relevant [po|j37[ . Equations with pairs of opposite families include 
sine-Gordon, mKdV and Broer-Kaup (two families each), NLS (four families). 



5.3 The singular manifold method as a Darboux transformation 

As qualitatively described in Section @, the singular manifold method 
looks very much like a resummation of the Laurent series, just like the geo- 
metric series 

-t-oo 

^x^x^O, (93) 

becomes a finite sum in the resummation variable X = x/{l — x) 

1 
Y,X-\X^O. (94) 

j=o 



The principle of the method is the following |119|. One first notices that 
the (infinite) Laurent series ( pq ) in the variable (p ~ (po can be rewritten as 
the sum of two terms 

u = VLogT + regular part. (95) 
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The first term PLogr, built from tlie singular part operator defined in Sec- 



tion 4.2, is a finite Laurent series and, if r is any variable fulfilling the two 



requirements for an expansion variable enunciated in Section 4.1, it captures 
all the singularities of u when ip —^ ipQ. The second term, temporarily called 
"regular part" for this reason, is yet unspecified. The sum of these two terms 
is therefore a finite Laurent series (hence the name truncated series), and the 
variable t is a resummation variable which has made the former infinite series 
mif — ifo'd, finite one. One then tries to identify this resummation (pB) with 
the definition of a Darboux transformation ([46|). This involves two features. 
The first feature is to uncover a link ([47] ) between r and a scalar component 
^ of a Lax pair. The second feature is to prove that the left over "regular 
part" is indeed a second solution to the PDE imder study. 



5.4 The degenerate case of linearizable equations 

The Burgers equation ([7l|), under the transformation of Forsyth (Ref. [Q 
p. 106), 

ti^aLogr, T — tp, (96) 

is linearized into the heat equation 

b^t + ^.cc + Git)^P = 0. (97) 



This can be considered as a degenerate Darboux transformation ([4 61), in 
which U is identically zero and ip satisfies a single linear equation, not a pair 
of linear equations, so this fits the general scheme. 

Another classical example is the second order particular Monge- Ampere 
equation s + pq ^ 0, linearized into the d'Alembert equation s = : 

s +pq = Uxt + UxUt = 0, (98) 

u = LogT, T = tp, ip^t = 0. (99) 

5.5 Choices of Lax pairs and equivalent Riccati pseudopotentials 

To fix the ideas, we list here a few usual second order and third order Lax 
pairs depending on undetermined coefficients, together with the constraints 
imposed on these coefficients by the commutativity condition. 

It is sometimes appropriate to represent an n-th order Lax pair by the 
2(n — 1) equations satisfied by an equivalent {n — l)-component pseudopo- 
tential Y of Riccati type, the first component of which is chosen as 

Yi = ^Ji>, (100) 

in which ip is a, scalar component of the Lax pair. 
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Second-order Lax pairs and their privilege 

The general second-order scalar Lax pair reads, in the case of two inde- 
pendent variables (x,t), 

Lilp = Ipx^j; - dVa; - 0,4' == 0, 
L21P = ipt - b-tpx - ctp ^ 0, 

Xq = -at + a^b + 2abx + c^, 
Xi = -dt + (6^ + 2c - bd)^ 

For the inverse scattering method to apply, the coefficients {d, a) of the a;-part 





(101) 




(102) 




(103) 




(104) 


- Cxd = 0, 


(105) 


0. 


(106) 



(101) are required to depend linearly on the field U of the PDE. 



The Lax pair (101)- (102) is identical to a linearized version of the Riccati 
system satisfied by the most general expansion variable Y defined by (pq), 
under the correspondence 

r = S-i-f-, B^Q, (107) 

Vx 

d = 2A, a^A^~A^- 5/2, b^-C, c = C^/2 + AC + j Atdx, (108) 

and the commutator of the Lax pair is (p3). 

In particular, when the coefficient d is zero or when, by a linear change 

/ ddx/2 
%l) \-^ eJ V"! it can be set to zero without altering the linearity of a on 

U, the correspondence is Q 

X^^, B = l, A = Q, (109) 

Wx 

d = 0, a = -5/2, 6 = -C, c = C^/2, (110) 

Li4j = 4Jxx + ^i^ ^ Q, (111) 

LzV- = ^t + C4^x - ^V' = 0, (112) 

2[ii, L2] = X^St + C^xx + CS, + 2a,S = 0. (113) 

Therefore second order Lax pairs are privileged in the singularity ap- 
proach, in the sense that their coefficients can be identified with the elemen- 
tary homographic invariants 5, C of the invariant Painleve analysis and, if 
appropriate. A, B. Conversely, and this has historically been the reason of 



some errors described in Section |8.2| , at the stage of searching for the BT, 
these homographic invariants 5, C are useless when the Lax order is higher 
than two and they should not be considered. 
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As explained in Section 5.3, given a Lax pair, one should define from it 
either one or two scalars ipf. Consider the second order Lax pair defined by 
the gradient of Y . Then, for one-family PDEs, this unique scalar ip is defined 
by (107). For two- family PDEs, the two scalars -0/ are defined by 



Y 



V'2' 



which leads to the zero-curvature representation of the Lax pair 



(*-m^; 



{dt - M) 



0, L- 



= 0, 



-A- B-^B^/2 B-^ 

B{A^-A^-S/2) A + B-^B^/2 



(114) 

(115) 
(116) 



M = 



AC + C^/2- B-'^Bt/2 



'CB- 



B{{CS + C^^)/2 + At + CA^ + C^A) -AC - C^/2 + B~^Bt/2 



The reason why the Riccati form is the most suitable characterization of the 
Lax pair is that it allows two linearizations |g8,101|, namely (107) and (114), 
depending on whether the PDE has one family or two opposite families. 



Third-order Lelx pairs 

The general third-order scalar Lax pair is defined as 

Lilp = Ipxxx - fipxx - aTpa; - b-ip = 0, (117) 

L21P = tpt~ cipxx - dipx - eip = 0, (118) 

[Li,L2] = Xo + Xidx+X2dl (119) 
-'^o = -bt - acx + Cxxx + bxxC + 2bcfx + bcxf - e^xf 

+3bcx^ + 3b^c^ + 3bdx + b^d = 0, (120) 
Xi = -at + Scxx + 2bxC + GxxC + d^xx + Sac^x + 2adx 

-l-Sa^Cj, -I- ibc^ + a^d + 2acf.j, + ac^f - 2e,.f - d^^f = 0, (121) 

X2 = -ft + [cf + cf^ + 2c^f + df + 2ac + c^^ -|- M^ + 3e)^ = 0.(122) 

An equivalent two-component pseudopotential is the projective Riccati 
one Y = (Yi, Y2) | |8^ , p7[ (written below, for simplicity, in the case f — 0) 

Yi ^^,Y2 = %, (123) 

Yi.^ = -Y^ + Y2, (124) 

Y2.x = -YiY2 + aYi+b, (125) 

Yi,t = -{dYi + cY2)Yi + {ac + d,)Yi + (c, + d}Y2 + e, + 6c (126) 
^{cY2 + dYi+e)^, (127) 

Y2,t = -{dYi + cY2)Y2 + {2ac^ + a^c + bc + d^^ + ad + 2e^)Yi 
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+ {c^^ + 2d^ + ac)Y2 + 2bc^ + b^c + bd + e^^, 


(128) 


\,tx - Yi^xt = Xi + X2Y1, 


(129) 


2,t£c — Y2^xt = —Xq + X2y2- 


(130) 



When there is no reason to distinguish between x and i, for instance 
because the PDE is invariant under the permutation (Lorentz transformation) 

V: {d^,dt)^{dud^), (131) 

it is natural to consider the following third-order matrix Lax pair, invariant 



under ( |l3l| ), defined in the basis {'ipx,4't,^) |37 




/i /2 fa \ 
= 0, L = ( gi 52 53 , (132) 

10 0/ 

/ 51 52 53 \ 
= 0, M ^ [ hi h2 hs] . (133) 

\Q 1 I 

In the equivalent projective Riccati components (^1,^2) 

Yi ^^,¥2 = ^, (134) 

with the property Yi.t — 12,2:, it is defined as 

Yux = -Y^ + fiYi + f2Y2 + /3, (135) 

>2,x = -Y1Y2 + giYi + 52^2 + 53, (136) 

Yi,t = -Y1Y2 + giYi + 52^2 + 53, (137) 

Y2,t^-Yi + hiYi + h2Y2 + h:i. (138) 

The nine functions fj , gj,hj, j = 1,2,3, must satisfy six cross-derivative con- 
ditions Xj ~ 

{Yi,x)t - {Yi,t)x ^Xo+ XiYi + X2Y2 = 0, (139) 

{Y2,x)t ~ {Y2,t)x = ^3 + X4Y1 + X5Y2 = 0, (140) 

easy to write exphcitly. It is worth noticing that there exists no such invariant 
second-order matrix Lax pair. 

5.6 The admissible relations between t and ip 

By elimination of dt , one of the two PDEs defining the BT to be found can 



be made an ODE, e.g. (p4) or (152). This nonlinear ODE, with coeflficients 
depending on U and, in the l-|-l-dimensional case, on an arbitrary constant A, 
has the property |90|1 of being linearizable. This very strong property restricts 
the admissible choices ( |47|) to a finite number of possibilities, and full details 
can be found in Musette lecture |85| . 
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6 The algorithm of the singular manifold method 

We now have all the ingredients to give a general exposition of the method 
in the form of an algorithm. The present exposition largely follows the lines 
of Ref. 0. 

The various situations thus implemented are : one-family and two-opposite- 
family PDEs, second or higher order Lax pair, various allowed links between 
the two sets of functions (t, -0). 

Consider a PDE (BSh with only one family of movable singularities or 
exactly two families of movable singularities with opposite values of uq, and 
denote V the singular part operator of either the unique family or anyone of 
the two opposite families. 

First step. Assume a Darboux transformation defined as 

u = U + P(Log n - Log ra), E{u) = 0, (141) 

with u a solution of the PDE under consideration, U an unspecified field 
which most of the time will be found to be a second solution of the PDE, Tf 
the "entire" function (or more precisely ratio of entire functions) attached to 
each family /. For one- family PDEs, one denotes ti = r, r2 = 1, so the DT 



assumption (141) becomes 

u = U + VLogT, E{u) ^0. (142) 



A consequence of the assumption (141) is the existence of the involution 



V/: (u,f/,T/)^-.(t/,«,r-i), (143) 

since the operator V is linear, and, for two-family PDEs, of the involution 

y{u,U): (I?,Ti,r2)^(-P,T2,Ti). (144) 

Second step. Choose the order two, then three, then . . . , for the unknown 
Lax pair, and define one or two (as many as the number of families) scalars 
ijjf from the component(s) of its wave vector (e.g. the scalar wave vector if 
the PDE has one family and the pair is defined in scalar form). Such sample 



Lax pairs and scalars can be found in Section 5.5. 
Third step. Choose an explicit link F 

V/: VLogTf ^ Fiijf), (145) 

the same for each family /, between the functions r/ and the scalars ^pf 
defined from the Lax pair. According to Section |5.6|, at each scattering order. 



there exists only a finite number of choices (145), among them the most 
frequent one 

V/: Tf=^f. (146) 
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Fourth step. Define the "truncation" and solve it, that is to say : with the 



ass ump ti ons (|l4l| ) for a DT, ( |145| ) for a hnk between r/ and ■0/, ( |l0lD -( |l02D 
or (117)-(118) or other for the Lax pair in -0, express E(u) as a polynomial in 
the derivatives of ijjf which is irreducible modulo the Lax pair. For the above 
pairs and a one- family PDE, this amounts to eliminate any derivative of tp of 
order in {x,t) higher than or equal to (2,0) or (0,1) (second order case) or 
to (3, 0) or (0, 1) (third order), thus resulting in a polynomial of one variable 
tpx/'ip (second order) or two variables '>px/'4','>Pxx/'>P (third order) 

E{u) = Y^ E^{S, C, U)i^p/i:xY+'^ (one-family PDE, second order), (147) 

k>0 1>Q 

(one-family PDE, third order). (148) 

Since one has no more information on this polynomial E{u) except the fact 
that it must vanish, one requests that it identically vanishes, by solving the 
set of determining equations 

Vj Ej{S, C, [/) = (one-family PDE, second order) (149) 

Vfc yi Ekj {a,b,c,d,e,U) = (one-family PDE, third order) (150) 

for the unknown coefficients {S, C) or (a, 5, c, d, e) as functions of U , and one 
establishes the constraint(s) on U by eliminating {S,C) or {a,b,c,d,e). The 
strategy of resolution is developed in Section |7.3| . 

The constraints on U reflect the integrability level of the PDE. If the only 
constraint on U is to satisfy some PDE, one is on the way to an auto-BT if the 
PDE for U is the same as the PDE for u, or to a remarkable correspondence 
(hetero-BT) between the two PDEs. 

The second, third and fourth steps must be repeated until a success occurs. 
The process is successful if and only if all the following conditions are met 

1. U comes out with one constraint exactly, namely : to be a solution of 
some PDE, 

2. (if an auto-BT is desired) the PDE satisfied by U is identical to (p3|), 

3. the vanishing of the commutator [Li, L2\ is equivalent to the vanishing 
of the PDE satisfied by U, 

4. in the l-fl-dimensional case only and if the PDE satisfied by U is identical 
to (pq), the coefficients depend on an arbitrary constant A, the spectral 
or Backlund parameter. 

At this stage, one has obtained the DT and the Lax pair. 

Fifth step. Obtain the two equations for the BT by eliminating ipj [Q 
between the DT and the Lax pair. This sometimes uneasy operation when the 
order n of the Lax pair is too high may become elementary by considering 
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the equivalent Riccati representation of the Lax pair and eliminating the 
appropriate components of Y rather than ■(/;. Assume for instance that t — ip, 



V = dx, and the PDE has only one family. Then Eq. (141) reads 



Yi=u-U (151) 



with Yi defined in (100), and the BT is computed as follows : eliminate all 



the components of Y but Yi between the equations for the gradient of Y, 



then in the resulting equations substitute Yi as defined in (151) 



If the computation of the BT requires the eUmination of ^2 between (124) 



(128), this BT is 



Yi,xx + 3YiYi,x + Y^^~aYi~b^0, (152) 

Yi,t - (cYi^. + cY,^ + dYi + e)x = 0, (153) 

{Yi,xx)t - {Yi,t)xx = ^0 + XiYi + X2YI = 0, (154) 



in which Yi is replaced by an expression of u — C/, e.g. (151). 

Although, let us repeat it, the method equally applies to integrable as 
well as nonintegrable PDEs, examples are split according to that distinction, 
to help the reader to choose his/her field of interest. 

6.1 Where to truncate, and with which variable? 

This Section is self-contained, and mainly destinated to persons accus- 
tomed to perform the WTC truncation. Although some paragraphs might be 
redundant with Section O, it may help the reader by presenting a comple- 
mentary point of view. 

Let us assume in this Section that the unknown Lax pair is second order. 
Then the truncation defined in the fourth step of Section p^ is performed in the 



style of Weiss et al. [119 , i.e. with a single variable. This WTC truncation 
consists in forcing the series (p5]) to terminate; let us denote p and q the 
singularity orders of m and E{u), ~p' the rank at which the series for u stops, 
and ^q' the corresponding rank of the series for E 



u = Y^ UjZ^+P, uau^p, ^0, E = J2 EjZ^+'^, (155) 

in which the truncation variable Z chosen by WTC is Z = (p — tpQ. Since one 
has no more information on Z, the method of WTC is to require the separate 
satisfaction of each of the truncation equations 

Vj-O,...,-?' :E,=0. (156) 

In earlier presentations of the method, one had to prove by recurrence 
that, assuming that enough consecutive coefficients Uj vanish beyond j = — p'. 
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then all further coefhcients Uj would vanish. This painful task is useless if 
one defines the process as done above. 

The first question to be solved is : what are the admissible values of p' , 
i.e. those which respect the condition u-pi ^ 0? 

The answer depends on the choice of the truncation variable Z . In Section 



4.1, three choices were presented, Z = cither tp — ipQ, x ot Y, respectively 
defined by equations (JSJ), (pSJ), (|6q), with the property that any two of their 
inverse are linearly dependent. 

The advantage of x or y over if — tpQ is the following. The gradient of x 
(resp. Y) is a polynomial of degree two in x (resp. Y), so each derivation of a 
monomial aZ'^ increases the degree by one, while the gradient oi ip — ipo is a 
polynomial of degree zero in ip — ipo, so each derivation decreases the degree 
by one. Consequently, one finds two solutions and only two to the condition 



u-p. ^ jioq] 



1. p' — p,q' — q, in which case the three truncations are identical, since the 
three sets of equations Ej = are equivalent (the finite sum ^ Ej Z^^'^ 
is just the same polynomial of Z~^ written with three choices for its base 
variable), 

2. for X S'Hd ^ only, p' — 2p, q' — 1q^ in which case the two truncations are 
different since the two sets of equations Ej ~ are inequivalent (they are 
equivalent only HA — 0). 

To perform the first truncation p' ^ p,q' = q, one must then choose 
Z = X since Y brings no more information and ip— ipo creates equivalent but 
lengthier expressions. 

To perform the second truncation p' = 2p, q' = 2q, one must choose 
Z — Y , since x would create the a priori constraint A ^ 0. 

The second question to be solved is : given some PDE with such and such 
structure of singularities, and assuming that one of the above two truncations 
is relevant (which is a separate topic), which one should be selected? 

The answer lies in the two elementary identities [pl| 



tanhz — = — 2isech 

tanhz 



2^ + ^- 



tanh z H — = 2 tanh 

tanhz 



2^ + ^- 



(157) 

Let us explain why on two examples, the ODEs whose general solution is 
tanh(a; — xq) and sech(a; — xq), namely 

E = u' + u"^ - I = 0, u = tanh(x - xq), (158) 

E = v''^ + a''^v^ - ■i;2 ^ 0, v = a sech(x - xq), (159) 



(this is just for convenience that we do not set a = 1). Equation ( |15§| ) has 
the single family 

p = —1, q = —2, uq = l,Fuchs indices = (—1), (160) 
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and equation (159) has the two opposite families 

p — —1, q ~ —4, vq — ia,Fuchs indices = (—1), (161) 

in which ia denotes any square root of —a^. The first truncation 
-p -9 

« = E "^■^'^''' ^ = E ^j^'^'' ^^ ■ ^1 = 0' (162) 

generates the respective results 

u = x'\S^-2, (163) 

V = iax^^, E2 = a^il - 5) = 0, £^3 = 0, E^ = -a^S'^/A, (164) 



thus providing (after integration of the Riccati ODE (p4[)) the general solution 



of equation (158), and no solution at all for equation (|15' 




The second truncation 

-2p -2q 

u = Y^ UjY^+P, E = J2 EjY^+\ Vj : Ej = 0, (165) 

i=o i=o 

generates the respective results 

u = B^^Y^^ + (l/4)Br, ^ = 0, S ^ -1/2, B arbitrary, (166) 

V = iaB-^Y-^ - {l/4)iaBY, A = 0, S = -1/2, B arbitrary, (167) 



thus providing, thanks to the identities (157), the general solution for both 
equations. 

The conclusions from this exercise which can be generalized are : 

1. for PDEs with only one family, the second truncation brings no additional 
information as compared to the first one and is always useless; 

2. for PDEs with two opposite families (two opposite values of uq for a same 
value of p), the first truncation can never provide the general solution and 
can only provide particular solutions, while the second one may provide 
the general solution. 

This defines the guideline to be followed in the respective Sections and 
y. The question of the relevance of the parameter B, which seems useless in 
the above two examples, is addressed in Section bl 



7 The singular manifold method applied to one-family 
PDEs 
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7.1 Integrable equations with a second order Lelk pair 

There is only one truncation variable, which must be chosen as %. 
Weiss introduced a nice notion, initially for one-family integrable equa- 
tions with a second order Lax pair, later extended to two-family such equa- 



tions by Pickering |101]. This is the following. 



Definition 21. (|20|) Consider the set oi—q+p determining equations (149) 
Ej — 0, which depend on (5, C, U). One calls singular manifold equation 
(SME) the result of the elimination of U between them. 

In the two-family situation, these determining equations also depend on 



{A,B), see (165), and the extension of this definition [101| is to also require 
the elimination of {A,B). 

Despite its name, originally restricted to integrable equations, the SME 
can be made of several equations in the nonintegrable case. 

The SME has the following properties. 

1. unicity, whatever be the integrability of the PDE, 

2. invariance under homography by construction p9| , i.e. dependence only 
on one Schwarzian S and as many C quantities as independent variables 
other than the one in the Schwarzian, 

3. the SME set is made of one and only one equation if and only if the PDE 
is integrable. 

Although one can define a SME whatever be the order of the Lax pair. 



it is inconsistent, as will be explained in Section B.2, to do so whenever this 
order is higher than two. 

The Liouville equation 

It is convenient to consider, following Zhiber and Shabat |l25|] , the equa- 
tion 

E{u) = u^t + ae" -I- aie~" + aae"^" = 0, a ^ 0, (168) 

which has the advantage to include the Liouville equation ai = 02 = 0, the 
sine-Gordon equation (oi 7^ 0, 02 = 0) and the Tzitzeica equation (ai = 
0, a2 7^ 0). As to the case aia2 7^ 0, it fails the test. Let us consider here the 
Liouville case. The results to be found are its auto-BT |^ and its hctero- 
BT to the d'Alembert equation. This will be achieved with two different 
truncations. 

Although not algebraic in u, the PDE is algebraic in either e" or e~". 



Equation (168) always possesses the family 

e"--(2/a)(^,(^t(¥.-^o)"', indices (-1,2), P = (2/a)a,9t. (169) 
For Liouville, this is the only family. 
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The special form of Liouville equation allows the assumption 

e" = VLogT + e^, E{u) = 0, X> = {2la)d^dt, (170) 

to be integrated twice to yield 

2 

u = -2LogT + V, E{u) ^^EjT^-^ = 0, (171) 

in which nothing is assumed on V . 

The Liouville equation is nongeneric for the singular manifold method in 
the sense that it is linearizable into another equation (thus, it should even 
not be part of the Section ^J). 

Therefore we define the first truncation in an exceptional way, namely we 
do not assume any linear relations on r = ^ and just treat t as the truncation 
variable. The three determining equations are then quite simple p7[ | 

Eo = 2T,Tt + ae^ = 0, (172) 

El = T^t = 0, (173) 

E2 = 14t = 0, (174) 

and their general solution depends on two arbitrary functions of one variable 



r = /(x)+5(i), (175) 

--r.n = -^f{x)g'(t), 
a a 

2 T^n 2 r{x)g'{t) 



e^ = --r.n = --f{x)g'(t), (176) 

a a 



aU{x)+g{t)r 



(177) 



2 T^Tt 



e^ =T-^e^ +-^=0. (178) 



a T 



Thus, the two fields u and V are the general solution of, respectively, the 
Liouville and d'Alembert equations. The hetero-BT between these two equa- 



tions is provided by the e limin ation of / and g between ( 176 ), ( 177 ) and the 



X— and t— derivatives of (171) 



(u - w), = aAe("+")/^ (179) 

(u + «)J = _2A-le("-"'/^ (180) 

in which v is another solution of d'Alembert equation defined by 

e^ = (ArO-^e^ = --A-^47^- (181) 

a g'(t) 



Remar k. W hen performing the truncation ( 170 ), Tamizhmani and Lak- 
shmanan |109| already found e^ — 0, r^t = as a particular solution, while 
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the above truncation (171) proves it to be the general solution. Another dif- 
ference between the two truncations is the presence of a field V in (171), 
which allows us to find in addition the hetero-BT between the Liouville and 
d'Alembert equations. 

Let us now define the second truncation, by the assumption 



'2LogT + W, 



(182) 



and the link (146), with ijj solutio n of t he Lax pair (101)-(102). Introducing 
the Riccati variable Y defined by (107), this second truncation is equivalent 
tof§7l 



u = -2LogY + W, y-^ =B{x~^+A), 

4 

E{u) = ^i?,(5,C, Ai?,M^)i"-'"', Vj : E, = 0, 



(183) 



and its result is recovered from the truncation of sine-Gordon in Section 9.1 
by simply setting ai = 0. 



The AKNS equation 
The AKNS equation || 

E{u) = Uxxxt + '^a~^{2{ux - (i)uxt + [ut - l)uxx) = (184) 

admits the single fanfily 

p ~ —1, q = —5, uq = a, indices (—1, 1,4, 6), V — adx, (185) 



so the ass ump ti on fo r the DT is (142). Let us choose at second step the scalar 
Lax pair (|lll| )-(112) for ip, at third step the link (146) between r and t p. Then 
there are only three non identically zero determining equations ( |l49[) [ p3[ 

E2 = 4:aSC + 8{Ut - 7) - lGC{Ux - /3) = 0, (186) 

E3 = ~a{CSx + -iSCx) + leCxiUx - /?) " SUxt + ^CUxx = 0, (187) 
^5 = E{U) + {a/2){2SSt - CSSx - Sxxt - SxCxx) - 2Sx{Ut - 7) 

-AStiUx ~P)~ ^SUxt + 2{SC + CxxWxx = 0, (188) 



plus the ever present condition X — 0, Eqn. (113). Their detailed resolution 
for {Ux, Ut) is as follows. One eliminates Ut between E2 and E^ 



E-i 



E2,x = -iC{~AUxx + aSx) = 0, 



(189) 



discards the nongeneric solution C — 0,Ut = '-f,St = 0, introduces an arbi- 
trary function of t after one integration, and solves for Ux 



Ux-P=ia/A){S + 2X{t)). 



(190) 
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Then E2 is solved for Ut 

Ut-i^ aX{t)C. (191) 

The cross-derivative condition Uxt — Utx is solved for St 

St = 4A(t)C^ - 2A'(t). (192) 



Substituting Ux,Ut, St and Cxxx taken from (|62|) in E^, one obtains the con- 
dition 

E5 = 2aX{t)X'{t) = 0, (193) 

which introduces the spectral parameter as the arbitrary constant A. 
The solution for {Ux, Ut) is 

Ux~P= (a/4)(5 + 2A), C/t - 7 = aXC, (194) 

and the elimination of U defines the SME 

^-4A = 0. (195) 

L^x 

The solution for {S, C) is 

S = {^/a){Ux - /?) - 2A, C = [Ut - 7)/(aA), (196) 

and its cross-derivative condition 

X = E{U)/{aX) = (197) 

creates on the field U the only constraint that U satisfy the AKNS PDE. 
The BT is the resuh of the substitution x"^ = {u- U)/a in (|4|)-(|6|). 

The KdV equation 

The Korteweg-de Vries equation for u (^6|) is defined in conservative form, 
so it is cheaper to process the potential form 

E{v) = bvt + vxxx ~ {i/a)vl + F{t) = 0, v = Ux. (198) 

Its unique family is 

p = -1, q = -4, Wo = -2a, indices (-1, 1, 6), V = -2adx. (199) 

With the assumption 

v^V + VLogT, E{v) = 0, (200) 
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for the DT, the choice of the second-order scalar (111)-(112) for -0, the hnk 



(146) between r and ip, one generates the three determining equations 



E2 = -2a{bC + 2S)-12V^^0, (201) 

Es EE 2a{bC - S)^ ^ 0, (202) 

E^ = E{V) + 1^2 - ^^3,. - 0. (203) 

After one integration of £'3, the system {E2, E3) is solved for {S, C) 

S = -2A(f) - {2/a)V^, bC = 4A(t) - {2/a)V^, (204) 

in which X{t) is an arbitrary integration function. Then E4, as seen from 
its above written compacted expression, expresses that V satisfies the PDE. 
Last, the cross-derivative condition ( |62|) 

X = -2A'(t) - 2{E(y))^/{ab) = (205) 

introduces the spectral parameter as an arbitrary complex constant and 
proves that a Lax pair has been obtained for the conservative (not the po- 
tential) equation. This Lax pair can be written, at the reader's taste, either 



in the scalar representation (|lll|)-(112), with U = ¥3 



Li = dl- U/a - A, (206) 

L.2 = bdt + (4A - 2Ula)d,^ + U^^/a, (207) 

a[Li,L2] = bUt + U^.^ - {6/a)UU,, (208) 



or in the zero-curvature representation (115)-(|116|) 



^-(c/Aa J)' ^'"'^ 

M-b-'( -^-l"" 2C//a-4A\ 

^^-C/,Ja + 2(C//a + A)(C//a-2A) U^ja )' ^^'"^ 

or in the Riccati representation for w = x~^ (see (p4)-(p5[) and (c9)) 

uj^^-^--uj^={^ + \\-uj\ (211) 

^t = [,-Cuj + C,/2), = fe-i((2C//a - 4A)cj - U.,la)^. (212) 

This last representation is by far the best one, for it allows one to deduce 
immediately two quite important informations, namely the auto-Backlund 
transformation of KdV and the hetero-Backlund transformation between 



KdV and mKdV. Firstly, the substitution of the inverse relation of (200) 



w = (w - V)l{-2a) (213) 
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in (211)-(212) provides the auto-BT for the conservative form of KdV 



a{v + V)^:, = -2a^X + {v-Vf/2, (214) 

aibiv + V)t - 2F'it)) = -{v - V)iv ~ V)^^ + 2{V^ + v^V^ + vl), (215) 

after suitable differential consequences of the a;-part have been added to the 
i-part in order to suppress A and cubic terms in (Ell 



Secon dly, t he elimination of U between (211)-( [2i2|) leads to the mKdV 



equation ( 404 ) for w, with the identification w = atj, v — X] since conv erse ly 
t he e limination of uj leads to the KdV equation for U , the system ( |21l[) - 
(^12|) also represents the hetero-BT between KdV and mKdV ([Q Eq. (5.16), 

As to the SME, it results from the elimination of V between (£'2, E3, E4) 

&C-5-6A = 0. (216) 

Most of these results for KdV were found in the original paper of WTC |ll9t . 



Remark. The transformation (211) between w = auj and U is often called 



a Miura transformation, but it is really just one half of the hetero-BT. The 
advantage of the hetero-BT it that it is invertible, while the Miura transfor- 
mation as defined in the previous sentence is not. 

7.2 Integrable equations with a third order Lelx pair 

Let us process a few PDEs which possess a third order Lax pair, and let 
us first perform their one-family truncation with the (wrong) assumption of 
a second order Lax pair, because this often provides interesting results. 

The Boussinesq equation 

The Bou ssine sq equation (Bq) is often defined in a two-component evo- 



lution form [|124 | 



sBnr,, r-^ - / ut-r.^^0, (a, /3, e) constant, 

Let us consider its one-component "potential" form 

pBq(f ) = vtt + e^ {{vx + af + [13^ /2,)v^xx) ,^ = 0, u^v^, r = vt. (218) 
Equation ( pi8[ ) has only one family of movable singularities 

p=-l, q = -5, indices (-1, 1,4, 6), V = 2f3^d^, (219) 



and it passes the Painleve test [117|. Since ( |218| ) is a conservation law, the 
computations can be reduced by considering the "second potential Bq" equa- 
tion 

ppBq(w) = wtt + e^ {{wxx + of + (/3^/3)wj;j;j;j;) =0, u = V .^ = Wxx, (220) 
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whose single family is of the logarithmic type w ~ 2/3^ Log x 

p = 0", q = -4, indices(-l, 0, 1, 6), V = 2/3^ (221) 

Let us assume for the would-be DT the relation 

w = 2/3^ Logr + W, ppBq(u;) = 0, (222) 



and for the link between r and ■0 the identity ( |146| ) 



Let us first assume that ip satisfies the second-order scalar Lax pair ( 111 ) 



(112). This is equivalent to the usual WTC truncation in the invariant for- 



malism 

4 

ppBq(w;) = ^ii;,X^~" = 0, (223) 

3=0 

and this generates the three determining equations 

E2 = {'i/5)(3^e^S - 2C^ ~ 4,6^ {W^x + a) = 0, (224) 

E3 = ~2{Ct - CC^ - {/3hy3)S^) = 0, (225) 

^4 = (5^2 - S3,x)/2 + Cj + P^ppBq{W) = 0. (226) 

From the last equation E4 = 0, the desired solution ppBq(Vl^) = cannot 
be generic, so this second-order assumption fails to provide the auto-BT. 
However, it does provide another information, namely a hetcro-BT between 
the Boussinesq PDE and another PDE. Indeed, under the natural parametric 
representation of E3 (which, by the way, would be the SME if the second order 
were the correct one), 

S = 3zt- i{l3efzl/2, C = {Pefz^, (227) 

the field z, by the cross-derivative condition (|62[), satisfies the modified Boussi- 
nesq equation [ pl[ 

MBq(z) = ztt + {{pef li)z,.,^., + 2{l3efztz,., - 2{(3efzlz,, = 0. (228) 



Just like for the KdV equation (Section 7.1), this leads, after a short 
computation left to the reader, to the hetero-BT between the Boussinesq 
and the modified Boussinesq equations. 



Going to third order, the assumption (E23) and (14C), with ip solution of 



the scalar Lax pair (117)-(hia), generates 



ppBq(u;) = J2J2 Ek,iY^Yl, k + l<2. (229) 



fc=o 1=0 



These six determining equations Ek.i = 0, plus the three cross-derivative 
conditions X, = 0, j = 0, 1, 2, are solved as follows in the Gel'fand-Dikii case 
/ = 0: |3,|0| 
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£^02 = (/3e)2 - c2 = ^ c = I3e, 

Eii=d^O ^d = 0, 

E20 = 3(y^ +a) + 2l3^a = ^ a = -3(T4 + a)/{2(3^), 

Eio = eVxx - Pe^^Q ^ e^ ^ f3''^eVxx, 

Xi = 3Vxt + S^eVxxx + ^P^ebx =0^6 = g{t) - 3{p-^Vxx + p-^e-^Vt)/4, 

Xo = (3/(4£/32))pBq(y) =0 ^V satisfies the PDE (218), 

Eoo = 2p^g'{t) =0 ^ git) = A, 

(230) 
in wliicli A is an arbitrary constant. Tlie coefficients a, b, c, d, e are 

a =. -i3/2)p-\Vx + a), b = X- {'i/A)p-^Vxx - {■i/A)p-^e-^Vu 
c^(3e, d = 0, e^ (3-h{Vx+a), (231) 

Xo = (3/(4e/32))pBq(F), Xi = 0, X2 = 0, (232) 

and t hey defin e a third-order Lax pair of the potential Boussinesq equation 
(111) [|l23yi24 |8l. 



The BT is just (152)"(153) or equivalently, after substitution of Yi = 

(« - y)/(2/32), 

[v - V)xx + il3-^e-\v + V)t + 3/?-2(„ _ v){{v + V)^ + 2a) 
+p-*{v-Vf -?>(3^\ = Q, (233) 

(w + V)xx - p-^e-^v - V)t + p-^{v ~V)(v- V)x = 0. (234) 

The Hirota-Satsuma equation 

Defined as l^^l 

HS(w) = [wxxt + {&/a)wxWt]x = 0, a 7^ 0, (235) 

it is better processed on its potential form 

pHS(u;) = Wxxt + {&la)wxwt + F{t) =0, a ^ 0. (236) 



The second order assumption ( |lll[ )-( |ll2| ) generates the three determining 
equations 

E2 = -2aSC - 6Wt + QCWx = 0, 

£3 EE aSt + 2aSCx - &CxWx = 0, (237) 

£4 = pHS(W^) - aiS'^C + Sxt/2 + SCxx) - iSWt + 3(5C + C:,,)VK, = 0. 

In the generic case Cx 7^ 0, their general solution is unknown, in particular 
we have not succeeded to perform the elimination of (5*, C) to find the con- 
straint(s) satisfied by W . It is easy to eliminate W but this gives rise to two 
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equations for {S, C) 



St 



6VF, = 2a5 + a-f, 6Wt 



CSt 

'7l 



M23 

M4E 



0, 



(238) 
(239) 



CSt \ _ / S^ 

Cx J X \ dx / t 

1 - 6a-'F{t) - C-\ACSSt + CS^^t + 2CSt) 
- C^^{2CS^ + C^StSx - 2CCxxSxt) - 2C-^CCl,St = 0, (240) 



and their possible functional de pen dence is unsettled. Anyhow, the field W 
cannot be a second solution of d23£ 



The third order assumption ( |ll7| )-( |ri8[) , with the link (f46) and the trun- 
cated expansion 



w 



W + adx Log r, 



(241) 



generates seven determining equations ( |150| ). The y are e asily solved [g6| and 
their uniqu e sol ution defines the Lax pair ( ^68|) -( 269 ), with W a second 
solution of (236). 



The Tzitzeica equation 

The equation is defined by (168), i n th e case (oi =0,02 7^ 0). It posseses 
two families, the first one defined by (169), the second one by 



^/{l/a2)'Px'Pt{'p- (Po) \ indices (-1,2). 



(242) 



These two families are not opposite, but the second family is irrelevant be- 
cause the Tzitzeica equation has a one-to-one correspondence p9[ with a 
one-family equation, namely the potential form (236) of the Hirota-Satsuma 
PDE in the particular case F{t) = 0. This correspondence is obtained by the 
elimination of 02 in equation (|16^) 



F{t) = 0, e" = —wt 
aa 



e-^"(e^"Tzi(u)), 



/pHSH 

V Wt 



(243) 



The irrelevance of the second family is confirmed by the negative result of 
Weiss |11S,1^ obtained when performing a truncation on e~". 

All the truncations will accordingly take the same form (|l7ClD as for the 



Liouville equation, which implies that t is an object invariant under the 



permutation (131). Depending on the Lax pair assumption, the link between 
r and ip will be e ither the identity (case of a scalar ip invariant under the 
permutation (131)) or not (if the scalar -0 is not invariant, e.g. because the 
Lax pair itself is not invariant), as detailed below. 

Let us first assume a second order Lax pair. To the author's knowledge, 
one cannot define a scalar ^, linked to such a Lax pair, which, like t, would 
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be invariant under (131). This is probably the reason why the assumption 

)-(112) generates so 



111 



T = ip with ip solution of the noninvariant Lax pair (] 
intricate determining equations that their general solution has not yet been 
obtained pm: these equations are however consistent in the sense that one 
easily finds the particular exact solution 



a2a 



= 2cp{x- ct- xi,g2,A+ -^-3-) - 2cp(a; + ct- 0:2,52,^ 



"8^ 



),(244) 



depending on five arbitrary constants (xi,X2,c,g2,A) and representing the 
superposition of two traveling waves of opposite velocities. 

From this second-order WTC truncation, and with appropriate assump- 
tions, one can also find a particular solution which represents a binary Dar- 



boux transformation |108 



Let us now turn to the third order assumption. One can postulate either 
a Lax pair invariant under ( 131 ), such as t he m at rix p air ( |l32|) -(133), or a 
noninvariant Lax pair such as the scalar pair (117)-( |ll8| ). In the first case, one 
must assume the identity link t = ip, while in the second case the assumed 
link must be noninvariant. Both assumptions lead to a success |Q. Let us 
detail here the invariant assumption, i.e. a priori the simpler one. 



The truncation is defined by (170), the link (146), and the matrix Lax 
pair (|13|)-(|13|) 



3 3-fc 



E{u) = J2J2 ^fe' (/j- ' 9j ' hj , U)Y,^Y^ ,yk,l : Eki^O, (245) 



fc=o 1=0 



in which (Y i, I2 ) are the two components of the projective Riccati pseu- 
dopotential (135)-( |l38|) equivalent to the Lax pair. To these ten determining 
equations in U and the nine u nkno w n coe fficients, one must add the six cross- 
derivative conditions Xj = (139)-( |l4(]|) . 

During their resolution, one first proves that the product f2hi cannot 
vanish (otherwise 02 would be zero). This makes the sixteen equations alge- 
braically independent and equivalent to the fifteen differential relations 



fj,t,g3.x,g].t,h.j,^,gj,^t = P{{fk,gk,hk},k= 1,2,3), j = 1,2,3, (246) 

with P polynomials whose coefficients depend on U, Ux,Ut, Uxt, plus the sin- 
gle algebraic relation 



Eon = 02 - — (53 + 9192 + (a/2)e^)^ = 0. 



(247) 
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They are solved successively as [equations are referenced as in (246)-(247)] 



93,xt — {93,x)t 

9i,x — g2,t 

92,t 



E, 



■00 



EiU) - 0, 

3 go{x, t) : gi^ go,*, 52 = 5o,x, 

33 = -ae^ - go,x9o,t ~ 9a.xt, 



^W~^f^\ 



notation W = e^ + {2/a)go,xt, 



92,x 

9s,x 

f2,t 

hi,x 
9i,t 

93,t 
/l2,x 



/s = -y/a^W Vo5o,t - fi9o,x - gl.^ + go.xx, 
/i = W^/W + 2go,^, 

h2 = Wt/W + 2go',t-fo.t/fo, 

fo.x = 0, 

hs = 5o,t(/o.t//o - Wt/W - go,t) + go,u - ^/a^W-^go,x/fo, 

fo,t = 0, 

go.xt = 0. 

(248) 
The irrelevant arbitrary function go reflects the freedom in the definition 
(170) of T and can be absorbed by redefining r as Te^^° . Thus the solution 
is unique : the field U must satisfy the Tzitzeica PDE, and /o is an arbitrary 
nonzero complex constant A. Accordingly, one has obtained a Lax pair and 
a Darboux transformation. The equivalent projective Riccati representation 
of the matrix Lax pair is 



-Y^ + L/.Yi 



u 



Yi,x 

Y2,x = -Y1Y2 - ae 

Yi^t = -Y1Y2 - ae^ 

r2 



^Xe-"Y2 



(249) 
(250) 
(251) 
(252) 



Y2^t = -Yi + UtY2 + V^A-^e-^Fi, 
with cross-derivative conditions proportional to the Tzitzeica equation 

(yi,x)t - {Yr^t)x = YrE{V)., (r2,x)t - {Y2^t)x - Y2E{V). (253) 

This Lax pair is the rewriting in matrix form of the scalar triplet given by 
Tzitzeica HI^ 

(254) 
(255) 
(256) 



-Txx + VxTx + ^fa^Xe ^Tt = 0, 
-Ttt + Utn + ^X'^e-^Tx = 0, 
-Txt - ae^T = 0. 

The Lax pair admits by construction the involution |111,50[ 



(r,e^,A) 



2 T^Tt 



,-A 



(257) 



equivalent to 



(T,e^,A)^(l/r,e^ + I?LogT,-A), 



(258) 
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which defines another, equivalent, writing of the Darboux transformation 

e" = -e^ - -^. (259) 



Remark. Knowing these resuhs, one can also write this DT |122,^ as 



a difference of the two fields u — U m terms of the two components of a 
projective Riccati pseudopotential 

u^U + Log(-2A2yiy2 - 1), % = a-'^^X-'e-^^^Y,, (260) 



in a quite similar manner to the DT of Liouville and sine-Gordon (392). 
However, the field u is multivalued. 

In order to find the BT, one must now eliminate one of the two equivalent 
projective components, and this defines two possible, different, eliminations. 



In the first elimination, one takes I2 from (E49[) and substitutes it into 
the three remaining equations, which results in 

Y2 - (Yi,. + Y^ - U^Yi)e^/{^\), (261) 

ODE = Fi,,, + 3^1 Fi,, + y3 „ e-^{e^)^^Yi + a^X = 0, (262) 
PDE = Yi^t + e^ ((FiFi,. + Y^) - Y^^U^) /(V^A) + ae^ = 0, (263) 

(252) = -YiE{U) — i-ODE + (2Yi - U^ + a:,)PDE = 0, (264) 

y/a2X 

[ODE,PDE] = (Yi,,,)t - {Yi^t)r. = Yi{e^^E{U))^. (265) 



Only two of them are functionally independent, as shown by relation (264), 



but the commutator (265) of equations (p63)-(263) shows that this elimina- 
tion fails to generate the auto-BT of Tzitzeica equation. 

However, it does provide another result, which we now derive. The ODE 



(262) belongs to the classification of Gambler - this is the number 5, see 
Section |5.6| -, it is linearizable by the transformation Yi = dx Log ip into a 
third-order linear ODE, with the relation r — -0 b etween the two functions. 
This transformation also linearizes the PDE (|263| ), and the resulting linear 
system 

Txxx - (Uxx + UI)tx + y^aAr = 0, (266) 

-^/wiXTt + e^Txx - Uxe^Tx = 0, (267) 

which cannot be a scalar Lax pair of the Tzitzeica equation, is, in fact, the 



scalar Lax pair of the Hirota-Satsuma equation (|235|), see Section 7.2 



Txxx - {Q/a)wxTx -f ylr = 0, (268) 

Art - {2/a)wtTxx + {2/a)wxtTx = 0, (269) 



under the change of variables ( 243 ) 
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In the second elimination, one takes Yi from (250) and substitutes it into 
the three remaining equations 



Yi 



-iY2^. + ae^)/Y2, 



(270) 



ODE = FzFa.xx - 2Y^'^^ - ([/.Fs + 3ae^)y2,. 

+y^Xe-^Yi ~ o?e^^ = 0, (271) 

PDE = Y2Y2,t + Yi - UtY^ + V^A-i(a + e-^Fa.x) = 0, (272) 
(251) = E{U) + {d^ - ae^r2~^)PDE 



-VaJA^^e^^Fa^^oDE = 0, 
[ODE,PDE] = (^2,..)* - {Y2,th. 

- (3ae^ + U^Y2 + 3^2,. - r2a,)£;(C/). 



(273) 



(274) 



Only two of them are functionally indepen dent , as shown b y the r elat ion 
(273), and the vanishing of the commutator (274) of equations (271)- (|272| ) is 
equivalent to the vanishing of the Tzitzeica equation for U . This elimination 
therefore generates the auto-BT of Tzitzeica equation, by the substitution 



Y2 



(a/2) / (e" - e^) dx 



(275) 



into (P71D-(P72D 



The ODE part (271) of the BT is equivalently written as [Bll 



W^, 



Wx 



Wx 



2 7T^ + aJa2\- 



w~W 



2wxWx 



(276) 



The nonlinear ODE ( [271[ ) again belongs to the equivalence class of the 
fifth Gambler equation (G5), see section d.C, and its linearization 



with the notation Y2 = (a/2)(w - W),e^ ^ Wx,e^ = W^ 



Y, 



-a-^e-'^dxLogie^^) 



(277) 



transforms the two equations (271)-(p73) into the third-order scalar Lax pair 
of the Gel'fand and Dikii type (i.e. / = in (pl7|)-(pl8|)) 



CiP = ip.xx + {2Uxx - UDiJx + {{2Uxx - Ul)x/2 + V^aA)V' = 0, (278) 

(279) 



(280) 



[£, M] = iEdl + {2{e^)xE + Ex)dx 

H<^^)xEx + [iUxx - Ul)E, E - E{U). 



Thus, the noninvariant (under (131)) link between r and •0 that o ne w ould 
have had to postulate if one had chosen the scalar Lax pair ( 117 )-(llg) is a 
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posteriori provided by the linearizing formula (277) and the Riccati equation 



(251), this is the invertible transformation 

e^r = (e^7/;)„ e^V = -a~Vt, (281) 



and it clearly breaks the invariance under (131) 



The Sawada-Kotera and Kaup-Kupershmidt equations 

Because of their duality p5| , |115[ , it is convenient to introduce simultane- 
ously the Sawada-Kotera equation (SK) and the Kaup-Kupershmidt equation 
(KK). These are defined as 

SK(u) = Put+ [ u^xxx + —uu^.x + -^v? ) = 0, (282) 

pSK(w) = /Jwt -f v,j,,j,,j,,j,,j, ^ v^v^^^ -I- -^v1 = 0, (283) 



KK(u) = /3-Uf + u^^^^ ^ uu^^ -t- —ul -f — M^ ^ 0, (284) 

y ol aol ol J ^ 

pKK(w) = (3vt + Vxxxxx H VxVa:xx + TT^lx H — 7^x = 0: (285) 

a la a^ 

in which u denotes the conservative field and v the potential one, with u — v^- 
Both equations have the Painleve property |115| . Each of them has two 
families ^^ 

pSK, Fl : p = -1, vq = a, indices - 1, 1, 2, 3, 10, (286) 

pSK, F2 : p = -1, vq = 2a, indices - 2,-1, 1, 5, 12, (287) 

pKK, ¥\:p^ -1, vo = a/2, indices - 1, 1, 3, 5, 7, (288) 

pKK, F2 : p = -1, wo = 4a, indices -7,-1, 1, 10, 12. (289) 

The singular part operator T) attached to a given family is 2? = vqBx- The two 
families have residues which are not opposite, but fortunately each potential 
equation possesses in its hierarchy a "minus-one" equation |120| | 

6 
pSK_i : v^xt + —VxVt = 0, (290) 

a 

pKK_i : VtVrxt - -;vlt + -VxVt = 0, (291) 

4 a 

which has only one family (the first one is nothing else than the Hirota- 



Satsuma PDF, already processed in Section 7^). The equations SK and KK 



are t here fore to be considered as possessing the single family Fl, Fqs. (286) 



and (288) 



Let us assume the one-family DT (EOOl) and, succesively, the second-order 



scalar Lax pair (111)-( |112D , then the third-order scalar one (117)-( 118 ) with 
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the Gel'fand-Dikii simplification / == 0. As to the hnk between r and -0, at 
second order this is the identity, while at third order it can be, as outlined in 



Section 5.6 and detailed in Musette lecture |85l, either the linearizing trans- 
formation of the fifth Gambler equation or that of the twenty-fifth Gambler 
equation. 

Therefore, at the fourth step of the singular manifold method, for each 
PDE, one has only three possibilities to examine : order two and Riccati, 
order three and (G5), order three and (G25). This is done in the next two 
Sections. 



The Sawada-Kotera equation 



First truncation (order two and Riccati). The one-family truncation (147) 
with T ~ tp generates the three equations ^^ 

Ei = pC~ AS^ + %S^^ + mSV^/a ~ 180(14/a)2 - 30V;,,/a = 0,(292) 
E5 = -I3C^ - 2SS^ + S^^., + 30^,T4/a = 0, (293) 

Ee = pSK(y) + {SEi - S5,,)/2 + dS^iW^Ja - S^/2) = 0. (294) 

These equations possess two solutions [pOl , a nongeneric one S^ = 

S^-k^/2, C==c + co, c = fcV/3 + 2co/3, 

V/a = a^ - (c " co)t, fcVl2 + /3co/9, 33) - k^x/l2 

+ ((5(fc^ + /3co)fcV36 - 1253)*)//?, (295) 

in which C, is the Weierstrass function and {k, co, 33) are arbitrary constants, 
and a generic one 5*0; 7^ defined by the four equations 

T4 = a(/3a + 255, - 5,,,)/(305,), (296) 

Vt = ... (297) 

Ml = (G,/5,),, -G- SlGl/b + 25G,/5, = 0, (298) 

M2 = 29 terms = 0, also vanishing if G == 0, (299) 

in which G is defined by 

G = S^^, + 452 - /3G. (300) 

The general solution (5, G) of the system Mi = 0, M2 = has not yet been 
obtained, for the elimination of 5 or G is difficult. This difficulty refiects the 
fact that second order is not the correct order. Neverthel ess, these compli- 
cated equations admit the very simple particular solution |115| (it would be 
interesting to prove that this is the general solution), 

G = 0, T4 = a5/3, pKK(y) = 0, (301) 
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SO the field V in the DT assumption (200) satisfies a different PDE. namely 
the potential KK equation. This defines a hetero-BT between the conservative 
forms of SK and KK ^,^1,^ 



a{v + V/2)^ + {v~ Vf = 0, P{v + V/2)t 
pSK(w) = 0, pKK(F) = 



0, 



(302) 



(see Ref. |58| for the exact expressio n of t he t- part) . 

With G = 0, the hnear system (in)-(112) is a degenerate Lax pair for 
KK, since it la cks a spectral parameter. 

Still when (301) holds, the field x~^ satisfies a fifth order PDE, the Fordy- 
Gibbons equation, and the explicit writing of its hetero-BT with the SK 
equation is left to the reader. 

Second truncation (order three and (G5)). This assumption creates no a 
priori constraint on the coefficients (a, 6) of the spectral problem (117), and 
the linearizing transformation of (G5) is just the identity t — ^. This gen- 



erates six determining equations (150). The process is successful 



115 



lia22 



and V is found to be a second solution of pSK (notation C/ = 14 as usual) 



b 
a 



A, 
-6t//a 



6—dx 
a 



(303) 
(304) 

(305) 



i2 = I3dt 



18 



C/x 



9A a 



36- 



.U' 



.u. 



a:r-36A-,(306) 
a 



[Li,L2] = ep-^a-'^SKiU). 



(307) 



This is the Lax pair given by Satsuma and Kaup [ 106 



The BT results from the elimination of I2, which provides Eqs. (152) 
( p3|) for Yi = Y, 



Y^^ + 3YY, + r^ + 6(C//a)r - A = 0, 
/3rt-9[(A-2f/,/a)(y, + y2) 

+ i{XU/a - [U/afY) + {2/i){U^^/a)Y], = 0, 
l3{{Y^,)t - {Yt)..)lY = -(6/a)SK(t/), 

followed by the substitution Y = {v — V)/a, 

{v - V).,^/a + i{v - V){v + V)^/a' + {v - Vf/a"" 
(3{v ~ V)t/a - (3/2) [(w - V)^^^^/a 
+ {b{v - V){v + y),,, + Vo{v + V)^{v - V)^^)lo? 
+(15(i; - Vf[v - V\^ + 30(t; - V){v + Vf^lo? 
+30(i; - Vf(y + V),/^* ^ g(^ ^ Vf loi% = 0, 



(308) 
(309) 

(310) 
A = 0, (311) 



(312) 



a result due to Satsuma and Kaup |106| 
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The Kaup-Kupershmidt equation 

First truncation (order two and Riccati). The one-family truncation (147) 
with T — ip generates the three equations BQ] 

E2 = 15(S'/4 - 3X4/a) = 0, 

Ei = I3C/2 + 7^2/4 + 3S',./4 - 15(^14 - 

Eq = {S/2)Ei/a + {AP{CS + C) - S^ 



Vxxx)/oi - 

{2\/2)Sl - 



90(14/a) 



+ (15/4) (3^2 _ 2Sxx)Vx/a - ^bS{Vx/af + pKK(y) - 0. 



(313) 
= 0, (314) 

(315) 



As opposed to the (difficult) SK case, these equations are easy to solve and 
possess the unique solution |115| 

Vx = aS/l2, pSK(y) = 0, S.,^ + S^/A - /3C = 0. 



(316) 



This is a strong indication that the particular solution (301) of ( |296| )-( p99| ) 
should be the general one. One again recovers, by a nice duality, the hetero- 
BT between KK and SK. 

Second truncation (order three and (G5)). This generates thirteen deter- 
mining equations ( |l5C| ). This truncation fails and provides no solution at all 
(one determining equation is i?2,2 = 45/8 = 0), not even the one-soliton 
solution. Indeed, the one-soliton solution of Kaup corresponds to constant 
coefficients for the scalar Lax pair (117)-( [llq) w ith f — 0, and, with the 
above procedure, the only way to obtain it | pl , [l01 is to enforce the two ffist 
integrals Ki and K2 which result from the zero value of b, 

Ki = il^xx - aijj, K2 = fAx - «V^^ - '2i'>Pxx - a'ip)ip. (317) 

Third truncation (order three and (G25)). This assumption implies, among 
the coefficients (a, 6) of the spectral problem (117), the a priori constraint 

m 

b-ax/2 = X{t), (318) 

and the linearizing transformation of (G25) defines the link between r and -0 

Tx X{t) 



1P 



(319) 



r Yi,x + {l/2)Y^^~a/2' 

This generates fourteen determining equations (|l50| ) (i.e. the same order of 
magnitude as for the (G5) assumption) in the basis (V'^/V'j ''Pxx/i^—i{ipx/'4')^ + 
bxi^x/ibip) — a)/2); they are solved as follows {gk denotes an arbitrary inte- 
gration function, A an arbitrary integration constant) |M 



^1,5 
^0,5 

El A 

Xo 



c = 9A(t)//3, 

a — — 6T4/a, 

d={Wxxx/a-i6Vx/a)^)/P, 

X{t) = A independent of t, 

pKK(\/) = 0, 

e = g2{t) + {S6XVx/a + 72VxVxx/a^ + 3T4.../a)//3, 

b = 54(0 - 3T4a;/a + Vxxx/iSa)), 

9i = A, 



(320) 
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and the result is, with U = Vx, 

b^ X- 3Ux/a, (321) 

a = -6U/a, (322) 

^^^°^^ " ^W^-(l/2)(V./^)2 + 3(C//a)- ^^2^) 

Li=dl + 6-dx + 3— - A, (324) 



::i2 1 / o ^^2; , r)^ \ ci o ^ xxx 

a O? I a 



L2 = Pdt - 9\dl + 3-^ + 36^ ] d^ - 3 



-72^-36A-, (325) 

/3[Li, L2] = (6/a)KK(C/)a, + (3/a)KK(C/),. (326) 

This i s the Lax pair given by Kaup ||65| . The integration of the first-order 
ODE ( p23|) modulo the Lax pair yields the DT given by Levi and Ragnisco 

T - ^jJ^Pxx - {1/2)^1 + 3(C//a)V'', T, = X^p\ (327) 

Although the relation r^/ip'^ — constant is the same as in the case of KdV 
(see Eq. (4.14) in Ref. |119]), it cannot be taken as an a priori assumption, 
it is the result of the method. 

Starting from the vacuum solution U — 0, the general solution ip of Liip = 

f = 1, ^3 = A, (328) 



in which ci,C2,C3, K are arbitrary complex constants, leads by (141) to the 
one-soliton solution of Kaup |p3] 

u = ia/2)dl Log(2 + cosh(fc/2)(a; - {k/2)H/(3)), keU (329) 

for the choice (01,02,03) = {0,j^,—j),K^ — — fc^/12, which corresponds to 
the entire function 

T = -(fcVl2)(2 + cosh(fc/2)(x - (fc/2)V/9))e'''/^^^''+'''/^^'*^'^^ (330) 

Let us now obtain the auto-BT of KK, by an elimination. In order to 
perform this elimination easily, it is convenient to choose one of the two 
components of the pseudopotential Y so as to characterize the DT, 

KK:'("-^^=^=Z. (331) 

a T 
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The chosen equivalent system is the system satisfied by (li , Z) 



Yi 



—-■, ^ — 7 

%jj T 

Yi,x = -^1^2 + AZ-i - 3[//a, 

Zrr. = 21"! Z — Z , 



iiY^.t - [9AriV2 - (3C/,,/a + m{Ulaf)Y^ 
+iU^.,^/a + 72UU^/a^ + 9XU/a].,, 
pZt = [l8At//a + 9X^Z-^ + 9XY^^ 

+ (45(C//a)2 + 6{UrJa) - l8iUja)Yi 
+27(U/a)Y,' + (9/4)Y,^)Z]. 



9\^Z-^ 



(332) 

(333) 
(334) 

(335) 



(336) 



The BT then arises from the elimination of Y\ between ( |333| ), (334) and 
(pS™ (Eq. (335) must be discarded), which results in the two equations for 

Y^., - {■i/4.)Y^/Y + 3yr,/2 + rV4 + 6(C//a)r - 2A = 0, (337) 

(3Yt - (3/16) [3r^ + iWY^ + 3or2y,, + gy,,,, 
+30(^3 + 21;,) (y, + 4k/q) + 6or(y, + w^/af 

+30i;(i;, + 4T4,/a) + 20r(y,,, + AV^^^/a)\, = 0, 
P{{Yxx)t - {Yt)xx)lY = -(6/a)KK((7), 



(338) 
(339) 



followed by the substitution Y = 2{v-V)/a, ^ 

{v - VU/a - (3/4)(« - V)ll{a{v - V)) 

+ 3(w - V){v + V)^/a^ + {v- Vf/a^ - A = 0, 

P{v - V)tla - (3/2)[2(« - y),,,,/a + 60(« - Vf{v + V)^/a^ 

+ 12(« - Vf/a'' + (10(w -V){v + V)^^^ + 30(« + V)^{v - V), 
+ 15(« - V)^{v + F),,)/a2 + (30(t; - Vf{v ~ V)^^ 



+ 60(w - F)(u + F)^ + 15(w ~V){v- V)l)/a^ 



0. 



(340) 



(341) 



The simple form of the conservative equations ( 312 ) a nd (341 ) re sult s from 
the addition of suitable differential consequences of ( [31l| ) a nd ( p40| ). The 
a;-part of the B T ha s alrea dy be given by Rogers and Carillo [103 for A = 0. 
If we write (p3|)-(|334|) in the variables (Yi, Z2 = Z^^), 



Yi^x 

Z2,x 



-riV2 4 
-2riZ2 ■ 



AZ2 

hi, 



3C//a, 



(342) 
(343) 



both systems (124)-(125) and (p42)-(p43) are coupled Riccati systems, with 
the difference that the transformation from Yi to 'ipx/4' is a point transforma- 
tion whil e t he o ne from Z2 to 4'x/'>P is a contact one. Thus, the Riccati sys- 
tem (124)-(125) is in the classification of linearizable coupled Riccati systems 
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given by Lie (this is the projective one), while the Riccati system (342)-(343) 
is outside it. 

A minor open problem is to find a bilinear BT for SK equation (the one 
given in Ref. pO| contains a nonbilinear term). 

7.3 Nonintegrable equations, second scattering order 

Strictly speaking, nonintegrable equations have no associated scattering 
order. What is meant in the title of this section is that one assumes a given 
scattering order to process some nonintegrable PDEs. 

For algebraic PDEs in two variables, particular solutions in which {S, C) 
are constant are quite easy to find. They correspond to solutions u polynomial 
in tanh ^{x — ct — xq). The privilege of tanh is to be the general solution of 
the unique first order first degree nonlinear ODE with the PP, namely the 
Riccati equation tanh +tanh +S/2, in the particular case S = constant. 

A characteristic feature of nonintegrable equations is the absence of a BT. 
Therefore, the iteration of Section can only generate a finite number of new 
solutions pgj.pl, this will be seen on examples. 

The Kuramoto-Sivashinsky equation 

It is worth to handle this example in detail because it exhibits all the 
features of what should be done and more importantly of what should not 
be done when solving truncation equations. 

The equation of Kuramoto and Sivashinsky (||) (notation ji = 19/i') pos- 
sesses a single family [§J,|9| 

p — —3, q~—7, uq = 120t^, indices —1,6, ,1? = 60z^9^ + 60^'92.. 

(344) 
and the orthogonality condition at index 6 is satisfied. Since equation (0) is 
a conservation law, we therefore study it on its potential form 

E = vt + Y + l^'^xx + vvxxxx + G{t) = 0, u^Vx, (345) 

which has the unique family 

p — —2, q — —6, vq — — 60j^, indices — 1, 2, , V = 60j^i9| + 60/x'. 

(346) 
Although the no-log condition at index i — 2 is not satisfied, the -^-series 

V = —QQi>x~^ + 60//' Log-!/) + f 2 + 0(x), V2 arbitrary function, (347) 



in which the gradients of ij: and x are given by (|68|)-(|69|) and (|64|)-(|65D, 
contains one logarithm only, which cancels by derivation. 
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The one-family truncation assumption is 

vt = woX~^ + vix~^ + f02 Log V' + W2, vq2 constant (348) 



equivalent to a truncated series (—3 : 0) for u. Substituting ( |348D in ( p45D 
and eliminating any derivative of x and tp, one obtains 

6 

i=0 



Together with the identity (|62[ ), this defines a system of eight equations in 
the six unknowns (uq , f i , f 02 , 1^2 , S", C) . 

Equations j = 0, 1, 2 are s olve d for uq, wi , 1^02, ^2 exactly as in the Painleve 
test and yield the values in ( ^347[ ). The next five equations (j = 3, 4, 5, 6 and 
(H)) now read || 

E-i = l2Qu{^C + lbvSr,+V2,a:)^Q. (350) 

Ei = 60 {-Qv'^S^rx - ^v^S^ - 20fi'iyS + 2vC^ + ll^i'^) = 0, (351) 

E^ = 1^3 + 60 {^fi'C + 2Q^i'vS, - 2v^SS^ + i/^^^^^ 

-i^C,, + A2,x) = 0, (352) 

Ee, = E{v2) + 30 {^I'C., ~ \^[i!'^S - ^^'^^(205'^ + S^^) 

-v'^^S^ + -iSl + ^SSxx)) =0, 
X = St^C.j,^,^^2C.j,S^CS^ = Q. (353) 

The principles to be obeyed during the resolution are the following. 

1. Never increase the differential order of a given variable. On the contrary, 
solve for the higher derivatives in terms of the lower ones, and substi- 
tute the result, as well as its differential consequences, in the remaining 
equations. 

2. Never integrate a differential equation, unless it is just a total derivative. 
On the contrary, perform an algebraic resolution. 

3. Never solve for a function of, say, one variable as an expression in several 
variables. 

4. Close the solution by exhausting all Schwarz cross-derivative conditions. 

This computation is systematic, and its algorithmic version is known as 
the construction of a differential Groebner basis 1^,0 • 

The full system is split into [Ez-, E4,E^), independent of dt, and {Eq,X), 
explicitly depending on dt- The subsystem (E^jE^) is first solved, according 
to rule 1, as a Cramer system for {v2,x, Sxx)- After substitution of {v2,x, Sxx) 
and their derivatives in all the other equations, equation E^ is solved, accord- 
ing to rule 1, for Cxx and the result is recognized as being an x— derivative. 
This allows us to solve for Cx after the introduction of an arbitrary integra- 
tion function A of t only. As to {Eq,X), they are solved, according to rule 1, 
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as a Cramer system in variables involving only <-derivatives, namely (w2,t, St), 
for expressions independent of dt- To summarize this first stage, the original 
system is now equivalent to 

,.'2 

(354) 

(355) 
(356) 



(357) 



^XX 


-l''-\:''^>^'-'>- 


a 


-T«'+¥«+C<^^+^^'- 


V2,x 


= C-15iySa:, 


V2,t 
V 


h' 2i/^ j/^ 1/ 




+1WS' I C + 15C^, ^'^'^^Sl, 


St 


-SS<^+-'-C'--'-^-'-f-' 




-CS^ + hvSl. 



(358) 



One equation, and only one, namely (354), is an ODE. Integrating it as 



an elliptic ODE for S 29 1 would create useless subcases and complications 
and should, according to rule 2, not be done. This ODE should also not be 
replaced by its first integral, because the integrating factor Sx could be, and 
will indeed be, zero. According to rule 3, it is also forbidden to eliminate 
A(i) by solving e.g. (p55|) for it. The only thing to do is (rule 4) to close this 



solution by cross-differentiation. There are two such conditions: 

{Sxx)t - [St).. = -7^(3A2 + 308A + 5324) + §^(15A + 374)5 



-i^A' + ^^Sl + 2buSSl = 0, (359) 



{v2.tU - {V2..)t = -Ct - ^(3A + 22)C - ^-^SC + ^-^S^C 

+i^(3A + 71)5:, - 2bb^l'uSSx 
-IbQu^S^Sx = 0. (360) 

The latter is solved for Ct and provides a third cross-derivative condition 
{Cx)t - {Ct)x _ m" 



3 (BIA^ -t- 4028A -I- 47476) 

5//'2 
(lOlA + 2322)5 -I- ^(55A + 96)5^ 



1415^^ ^3 825^/^, 



52 R. Conte 

+ ^A' + ^-^Sl + 27bvSSl = 0. (361) 

This ends the hnear part of the resolution, and now conies the nonlin- 



ear part (algebraic Groebner). The two remaining equations (35£) and (361), 
considered as nonlinear in the two unknowns (Sx , S) , imply without compu- 
tation S't = 0, which then allows one to solve (t354) for the monomial S"^ as a 



polynomial in 5 of a smaller degree. Equations (|359| ) and (361) thus become 
linear in S and A', and their resultant in A' factorizes as a product of linear 
factors 

(A + 33)(A-ll + ^5) = 0. (362) 

/^ 

The first factor yields no solution. The second one provides the unique solu- 
tion 

S — -—— ] [ S + -—— 1=0, C = arbitrary constant c (363) 

38j^/ \ iv>v J 

and it leads to the two-parameter (c, xq) solution (0). The two equations 



( P63D represent the SME. 

If one performs the iteration of Section 0, starting from u = c, one gener- 
ates the solitary wave (0) and no more [E9| . 

The reduction u{x,t) = c + U{^),£, =^ x - ct oi the PDE (|) yields the 
ODE 

vU'" + nU' + C/V2 + K = Q, K arbitrary, (364) 



with the indices —1, (13 ± i\/7T)/2. Due to the two irrational indices, the 



general analytic solution (see definition in Section 9.3) can only depend on 
one arbitrary constant. This one-parameter solution, whose local expansion 
contains no logarithm, is known globally only for K = — 450j^A;^/(19^/i), this 
is (0), but its closed form expression for any K is still an open problem 



Being autonomous, the ODE (364) is equivalent to the nonautonomous 
second order ODE for V{U) 

an equation which has been studied from the Hamiltonian point of view H . 
7.4 Nonintegrable equations, third scattering order 



An example is given in the Section 9.3 



8 Two common errors in the one-family truncation 

Two errors are frequently made in the method of section [^. 
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8.1 The constant level term does not define a BT 

Consider the one-family truncation as done by WTC (the subscript T 
means "truncated") 

«WTC^^^WTC^,+p (3gg) 

i=o 

in which ip is the function defining the singularity manifold. 

In the WTC truncation, one considers three solutions of the PDE 

1. the Ihs u^'^'^ of the truncation p6q), 

2. the "constant level" coefficient u'^J-'^, 

3. the field U which appears in the Lax pair after the successful completion 
of the method. 

The frequently encountered argument "The constant level coefficient w^J"^ 
also satisfies the PDE, therefore one has obtained a BT" is wrong. This is 
obvious, since nonintegrable PDEs, which have no BT, nevertheless have this 
property. One can check it by taking the explicit example of a nonintegrable 
PDE |§. 

A hint that the above argument might be wrong is the fact, observed 
on all successful truncations, that the U in the Lax pair is never u^^'^. 
Let us prove this fact, with the homographically invariant analysis |l9|. The 
truncation of the same variable in the invariant formalism is 

UT = Y.u,x'^P, (367) 

3=0 

in which x is given by (p8|). This ut depends on the movable constant ipQ 
and one has 



WTC / 

u_p = UT((po = oo 



WTC _ „, /,„ \ (368) 



Since the results of the truncation do not depend on the movable constant (pQ , 
this proves that the Ihs u^"""*^ of the truncation and the constant level coeffi- 
cient u^p'~^ are not considered as distinct by the singular manifold method. 
Since the U in the Lax pair cannot be the truncated u (otherwise one would 
not have a Darboux transformation), this ends the proof. 

8.2 The WTC truncation is suitable iff' the Lax order is two 

We mean the truncation as originally introduced, not its updated version 
of Section 0. 

When the Lax pair has second order, everything is consistent. When the 
Lax pair has a higher order, e.g. three, the original method, as well as its 
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original invariant version p 
stage, it generates the —q 



, presents the foUowing inconsistency. In a first 
p equations Ej{S,C,U) 



of formula ( |147| ), 
which intrinsically correspond to a seconrf-order scattering problem (and this 
is precisely the inconsistency), and in a second stage it injects in each of 
these ~q + p equations a link between {S, C) and the scalar field ip of the 
Lax pair of higher order, thus generating determining equations which are 
hybrid between the second order and the higher one. The first nearly correct 
treatment has been made in Ref. K^. 

For the same reason, in order to obtain the Lax pair when its order is 
higher than two, it is also inconsistent to consider the so-called singular man- 



ifold equation (SME) |11£,19,101|, defined in Section 7.L When the Lax order 
is three, the correct extension of this SME notion would be the set of three 
relations on (a, b, c, d) resulting from the elimination of U between the four 
coefficients of the Lax pair (e is derivable from (122) so we discard it), but 
this seems of little interest. 

Although these inconsistencies may still provide the full result for some 



"robust" equations (Boussinesq |117|, Sawada-Kotera |115|, Hirota-Satsuma 
[p6|), there do exist equations for which it leads to a failure, and the Kaup- 
Kupershmidt equation |BG| is one of them. 



9 The singular manifold method applied to two-family 
PDEs 



By two-family, we mean two opposite families. 

This includes also the one-family truncation as a particular case. 

When the base member of the hierarchy of integrable equations has more 
than a single family, these families usually come by pairs of opposite singular 
part operators, just hke (P2)--(P6). Examples are enumerated at the end of 
Section 5.2, Then the sum of the two opposite singular parts 



V Log Ti—V Log T2 



only depends on the variable 



T2 



(369) 



(370) 



The current status of the method |8J,101|, which used to be called the 
two-singular manifold method |p3, is as follows. Most of the method for 
one- family equations still applies, with the difference that it is much more 
convenient to represent the Lax pairs in a Riccati fornr than in a scalar linear 
form. Let us restrict here to second-order scatte ring problems (for the third 
order case, see Section 9^) and to identity links (146) between the two t and 
the two ij: f unctions. Then Y satisfies a Riccati system and, as explained in 
Section 3.1, its most general expression is given by (p6|). 
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In the first step, r is simply replaced by Y in the assumption (141) for a 
DT. 

In the second step, the scattering problem is represented by the Riccati 
system satisfied by Y, whose coefficients depend on {S, C, A, B). 

The fourth step contains the main difference. Rather than truncating u 



at the level j = —p, one truncates it at the level j = —2p [B8 101], in order to 



implement the two movable singularities ti = and T2 = 0. So the truncation 



is [101 1 (for second order Lax pairs only) 



u^VLogY + U, (371) 

Y-'=B{x-'+A), (372) 

-2q 

i?(u)-^i?j(5,C,AS,C/)i"^+', (373) 

Vj Ej{S,C,A,B,U) = 0, (374) 

in which nothing is imposed on U. 

Let us remark that the relation A ^ does not characterize two-family 



PDEs, see the Liouville case in Section 9.1 



9.1 Integrable equations w^ith a second order Lelx pair 



The sine-Gordon equation 

The sine-Gordon equati on is defined for convenience as the case ai ^ 



0,02 =0 of the equation ( |168[) . Although not algebraic in u, it becomes 
algebraic in e" and it possesses two opposite families (opposite in the field 
u), both with p — —2, q = —2 

e" - -{2/a)ip^ipt{(p - fo)-^, indices (-1, 2), V = (2/a)d,dt. (375) 
e-"-(2/al)(/>^(^t(^-^o)~^ indices (-1,2), V ^ -{2/ai)d^dt. (376) 

The resulting DT assumption 

e" + (ai/a)e"" = {2/a)d:,dt LogF + W, E{u) = (377) 



with Y defined by ( p70D , can be integrated twice due to the special form of 
the PDE, resulting in 

u^-2hogY + W, E{u)^0, (378) 

in which nothing is imposed on W (we use W to reserve the symbol U for 
future use). For oi = 0, this truncation is what was called in Section ^^ the 
second truncation of Liouville equation. 
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The five determining equations in the unknowns (S*, C, A, B, W) are [101 3^] 

Eo = aB^e^ - 2C = 0, (379) 

El = 2{Ca, + 2AC) = 0, (380) 

E2 = 0, (Fuchs index) (381) 

^3 = ~<Tt ~ a{C., + 2AC) = 0, (382) 

E^ = a {Ca + (C^ + 2AC)^) /2 + a^B-'^e-'^ = 0, (383) 

with the abbreviation 

a = S + 2A^ -2A^, (384) 



and, together with the cross-derivative condition (pS), they are solved as 
usual by ascending values of j 



^0 


B^e"^ = -C, 
a 


(385) 


El 


A^-^{LogCU 


(386) 


E^ 


^--nx) + §-%, 


(387) 


Ei 


CC.,t - C.Ct + F{x)C^ + aiaF{x)-^C ^ 0, 


(388) 


X : 


aiF'{x) = 0. 


(389) 



in which i^ is a function of integration. For sine-Gordon, F{x) must be a 
constant 

F{x) = 2A2. (390) 

In the Liouville case, for which the truncation imposes no restriction on F{x), 
let us also require that F{x) be a constant. Then, for both equations, LogC 
is proportional to a second solution U of the PDE 

C=^X-^e^, EiU)^0, (391) 

and one has obtained the Darboux transformation 

u ^ ~2Logy + U, y = XBY, (392) 

in which y satisfies the Riccati system 

y:, = X + U.,y-Xy\ (393) 

2/* = -fA-i(e^-f(ai/a)e-V), (394) 

{Logy),t - (Log2/)t, = EiU). (395) 



Exact solutions by singularity analysis 19 February 2000 
The linearization 

yields the second-order matrix Lax pair 



[d. - L) 



■4^2 
^1 



0, L 



A 



'U^/2 



(a.-M)(j;)=0,M = -(a/2)A-(_(^^ 







-{ai/a)e 



e 
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(396) 

(397) 
(398) 



The auto-BT (classical for sine-Gord on, R e f. ||8l|] for Liouville) results 
from the substitution y = e-(«-f^)/2 into (^93|)-(|394|) 



[u + U)x = -4A sinh 
{u-U)t = \~^ (ae(" 






flie 



-(u+U)/2 



(399) 
(400) 



It coincides in the sine-Gordon case with the one given earlier, equations 
®-(0)- The ODE part ( p93| ) of the BT is a Riccati equation. 
The SME is |0l| 



S + C-^C - -C-^Cl + 2A2 ^ 0, 



(401) 



and it coincides, but this is not generic, with the one [116 19[| obtained from 
the (incorrect) truncation in x- 
Remarks. 

1. The reason for the presence of the apparently useless param eter B in the 
definition ( |66|) is to allow the precise correspondence (|l46| ) 



for some choice of i?, namely 

B = A"\ y = Y., W = U. 



(402) 



(403) 



2. In the Liouville case ai ~ 0, this is an example of a PDE with only one 
family and a nonzero value of A. 

The modified Korteweg-de Vries equation 

This PDE has the same scattering problem as sine-Gordon, so the com- 
putation should be, and indeed is, quite similar to that for sine-Gordon. 
Since this PDE has the conservative form 



iKdV(w) EE bwt + {w^^ -2{w- I3f/a^ + 6i/w)^ = 0, w = r^, (404) 
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it is technically cheaper to process its potential form 

p-mKdV(r) = brt + r,,, - 2(r, - /3)Va' + 6z/(r, - /?) + F{t) = 0. (405) 

Its invariance under the involution w — p <—> —{w — (3) provides an elegant 
way |7q] to derive the BT of the KdV equation and its hierarchy. Although 
the constants (3 and v could be set to zero by a transformation on (r, x, t) 
preserving the PP, it is convenient to keep them nonzero, for reasons explained 
at the end of this Section. This last PDE admits the two opposite families 
[a is any square root of a^) 

p~Q^, (7 = — 3, r^^aLogi/', indices (—1, 0,4), V — a. (406) 

The truncation is defined by 

r = aljOgY + R, (407) 



with ( |37l|)-( 374), and this generates five equations Ej = [101 , with the 



notation (384) for a 



Ei=QaA-Q{{R-ahogB)^-[3) = Q, (408) 

E2 = a{2A^ + AA^ -bC -2a + 6iy) - a'^Ei + 6aAf/6 = 0, (409) 
E3 = p-mKdV(i? - a Log S) - (3/2)a"V^ 

+ (cr - 4^2 - (l/3)a"^£;i^:r - '2.A^)Ex - 2AEi^^ - 2AE2 - E2.cc (410) 

E4 = expression vanishing with Ei,E2,E3,E^, (411) 

E5 = (3/4)aaa, + (1/4)^2^1 = 0, (412) 

X = = St + C,,. + 2CccS + CScc = 0. (413) 

They depend on {R, B) only through the combination R—a Log B. Equation 
j = A is a differential consequence of equations j — 1, 2, 3, 5, because 4 is a 
Fuchs index, and the other equations have been written so as to display how 
they are solved: 

Ei: A = a-\iR-aLogB),-l3), (414) 

E5 : a = -2(A(t)^ - i/), A arbitrary function, (415) 

E2: bC ^ 2A:, - 2A^ + 4:\{tf + 2z/, (416) 

E3 : p-mKdV(i? - a Log B) = 0, (417) 

X : X'(t) = 0. (418) 

Thus, their general solution can be expressed in terms of a second solution 



W of the mKdV equation (404) and an arbitrary complex constant A |101] 



W = iR^aLogB),c, A=iW~l3)/a, 
bC = 2W,,/a - 2{W ~ pf/a^ + 2v + 4A^ 
S = 2W^la-2{W -Pf/a^ + 2v~2X^, (419) 
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and the cross-derivative condition Xi = (Eq. (106)), equivalent to the 



mKdV equation (404) for W, proves that one has obtained a Darboux trans- 
formation and a Lax pair. 

The SME, obtained by the chmination of W between S and C, 

bC-S-6X^ = 0, (420) 



is identical to that of the KdV equation (216) 



The auto-BT of mKdV is obtained by the substitution 

Log{BY) = a-^ f{w - W)dx (421) 

in the two equations for the gradient of y = BY 

^^A(l-y)- 2^^1:1^, (422) 

y y oi 

5^^ir-4A^^l^ + (2(^-/)%2^-4A^),) . (423) 
y y \ a a^ a J ^ 

In the same manner as in the KdV truncation, these two Riccati equations 
can also be interpreted as the hetero-BT between the mKdV equation and the 
PDE satisfied by the pseudopotential y, called the Chen-Calogero-Degasperis- 
Fokas PDE. 

The nonlinear Schrodinger equation 

For the AKNS system of two second order equations in (u, v) (whose 
reduction u — v is NLS, see (|39|)), no two-family truncation has yet been 
defined which strictly follows the method and provides the desired result. 
It should be noted that the fourth order equation for u resulting from the 
elimination of w, known as the Broer-Kaup equation or classical Boussinesq 
system, admits a two-family truncation without any problem pq . 

The full result (DT, BT) can be found pi for the AKNS system by 



performing the one-family truncation |117] and then applying four involutions 
to the result of Weiss. 

A second open problem for this PDE is that its bilinear BT is not yet 
known. 



9.2 Integrable equations vifith a third order Lelx pair 

In principle, there is no additional difficulty to extend the method to a 
scattering order higher than two. A good equation to process would be the 
modified Boussinesq equation 

^^f-.* + (.-(3/2)aV) =0 

1 -vt - ia [Uxx ~ uv + a u )x ^0, 
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which has two opposite famiUes 



u^{2/a)x-\v^(ix~' (425) 

and a third order Lax pair hke the Boussinesq equation. 
The one-family assumption B^ 

u^U + {2/a)d^ Logr, v^V -Qdl Logr, E{u, v) = E{U, V) = 0, (426) 



with the identity hnk t — ip and the choice of the scalar Lax pair (117)-(118), 
already leads to the solution 

/ = -(3/2)aC/, a^{V - 3a^U^ - 3af/^)/4, 6 = A, 

c = -3a, d = -3a^U, e = 0. (427) 

Despite this success, it would be more consistent with the two-family 
structure to process this PDE with a two-family assumption, removing in 
passing the restriction E{U, V) ^ in ([426| ). This coul d make the coefficients 



(/, a, b) linear in ([/, V), which is not the case in (427). 

Table |l| summarizes, for a sample of PDEs, the currently best method to 
obtain its Lax pair, Darboux and Backlund transformations from a trunca- 
tion. 

9.3 Nonintegrable equations, second and third scattering order 

A nonintegrable equation has no determined scattering order, so this sec- 
tion cannot be split according to the scattering order. 

The KPP equation 

The KPP equation i jn\ ) possesses the two opposite families ([78|)-~(p0|) and 
it fails the test at index 4, so there can only exist particular solutions. Let us 
first review all the known solutions to this equation. 

In addition to the notation (|8l|), it is convenient to introduce the sym- 
metric constant 

ai = (2ei - 62 - e3)(2e2 - eg - ei)(2e3 - ei - ei)/{3df (428) 

and the entire function 

3d 



Ws^Y. ^-e''"^'' ^ (3/&)fc"*), fe„ = ^\-^, C„ arbitrary, (429) 



n=l 



i.e. the general solution of the third order linear system (117)-(118) with 
constant coefficients ^l| 

( fV^...-3a2^ -ai^ = 0, 

[b->pt- 3V'xx = 0. 
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Table 1. The relevant truncation for some 1 + 1-dimensional PDEs. The 
successive columns are : the usual name of the PDE (a p means the potential 
equation), its number of families (a * indicates that only one family is rele- 
vant, see details in Ref), the order of its Lax pair, the truncation variable(s), 
the link between r and tp, the singularity orders of u and E{u), the Fuchs 
indices (without the ever present —1), the number of determining equations, 
the reference to the place where the right method was first applied (earlier 
references may be found in it). The "?" in the AKNS system entry (the one 
whose NLS is a reduction) means that the method has not yet been applied 
to it, see text. 



Name 


f 


Lax 


Trunc. var. 


T 


-p ■■ -q 


indices 


nb. det. eq. 


Ref 


Liouville 


1 




r 




0:2 


2 


3 


3^] 


KdV 


1 


2 


X 


i, 


2 : 5 


4,6 


2 


119| 


AKNS eq. 


1 


2 


X 


i, 


1 : 5 


4,6 


3 


lol 


p-mKdV 


2 


2 


Y 


^ 


:3 


0,4 


4 


sine- Gordon 


2 


2 


Y 


^ 


: 2 


2 


4 


101 1 
101 1 


Broer-Kaup 


2 


2 


Y 


V 


:4 


0,3,4 


4 


pp-Boussinesq 


1 


3 


{i'./i>,i>../^) 


^ 


:4 


0,1,6 


6 


87 1 


p-SK 


1* 


3 


{■ipx/i>,i>xx/ip) 


^ 


1 :6 


1,2,3,10 


6 


9(| 


p-KK 


1* 


3 


(^./^,^../^) 


G25(i/;) 


1 :6 


1,3,5,7 


14 


Tzitzeica 


1* 


3 


{^b./^,^tm 


^ 


2:6 


2 


10 


37 1 


AKNS system 


4 


2 


7 




1 :3,1 : 3 


0,3,4 




3^1 



Let us also denote {j,l,m) any permutation of (1,2,3). Three distinct 
solutions are presently known. 

The first solution is trigonometric, this is a collision of two fronts |66|] 



U=J+dd., Logtf'g, C1C2C3 ^ 



(431) 



which depends on two arbitrary constants C1/C3, C2/C3. For Cj — 0, CiCm 7^ 
0, it degenerates into three heteroclinic (i.e. with different limits at both 
infinities) propagating fronts which depend on one arbitrary constant xq 



\- d- tanh -(X - ct- Xq), 



(432) 



k^ = [h - k^f, c = -3{ki + k^)/b. 

The second solution is elliptic |lO[ , 

u = si/3 -I- dipxy^p(tp), ^ = if'3, 53 = 0, 52 arbitrary, ai = 0, (433) 

it only exists under the constraint (codimcnsion is one) that one root Cj be at 
the middle of the two others and it depends on the four arbitrary constants 
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Ci,C2,C3,g2- Its degeneracy g2 = (i.e. p{ip) = ip~^) is the degeneracy 



ai = of the coUision of two fronts solution (431) 



The third and last solution is the stationary elliptic solution u(x) 

u{x) : -u" + 2d-'^{u-ei){u-e2){u-e3) = 0. (434) 

A trigonometric degeneracy bounded at infinity is made of the three homo- 
clinic stationary pulses 

u — Cj^ j= — sechi -j^{x — xq), ai — 0, 2ej — ei — Cm = 0,(435) 

V 2 dv 2 

it has codiniension one and it depends on the arbitrary constant xq. 

Let us now apply the various methods we have seen, in order to retrieve 
these solutions, namely 

1. enforcement of one of the two no- log conditions (|8^), 

2. enforcement of the two no- log conditions ([8^), 

3. one-family truncation with a second order assumption, 

4. one-family truncation with a third order assumption, 

5. two-family truncation with a second order assumption, 

6. two-family truncation with a third order assumption. 



The single no-log condition (82) has two solutions. The first one C = 



implies u* = and thus defines the reduction u{x), i.e. the elliptic equation 



(434). The second one is a first order nonlinear PDE for C{x,t), integrated 



by the method of characteristics as [^ 

F(/i, /2) = 0, c„ = 3fc„/6 = (3e„ - si)/(M), 

I, = e^C ~ c,r {C - C2T{C - c^Y\ (436) 

/2 = e*(C - ci)^i {C - C2T{C - C3)«^ 

in which p„,g„ depend on Cj. Unless some specific choice of the arbitrary 
function F is made, or (x, t) are no more taken as the independent variables, 
one cannot integrate further the system (|64|)-(|65|) for x and S. 

The two no-log conditions (p2| ) together, apart the already encountered 
solution C = 0, provide the two relations 

ax =0, ei = (62 -1-63)72, 

\?d^C'^ - (9/4)(62 - 63)^C - ihd^(Ct + CC^) = 0, (437) 

whose solution is similarly 

ai=0, F(/i,/2) = 0, 

_ 6C-(3(62-63)/(2d)) (e. - 63)^/d .438^ 

^ 6C+(3(62-63)/(2d))^ ' ^^"''^ 

^' =(62 - es)h/ibd') 



(foC)2~(3(62-e3)/(2d))2 
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The one- family truncation (142) with t — tjj and the second order as- 
sumption (111)-(112) generates three determining equations Ej{S,C,U) ~ 
0, j — 1, 2, 3. After solving the first one for U 

U = si/3 - bdC/6, (439) 

the two remaining equations, with the ever present condition (|62[), are pl| 

E2 = -602 -S + {b^/6)C^ - 6C^ = (440) 

^3 EE a2bC ~ 2ai + bSC/2 - b^C^ /108 

+S^/2-b^Ct/6 + 2bC^j3 = 0. (441) 

The elimination of S and Ct yields a factorizcd equation 

-b^Ct + bC^^ - 2b^CC^ + (4/9) (feC)^ - l2a2bC - 12ai = 0, (442) 
[&C:,:, - b^CC^c + 6^CV9 - 'ia2bC - 3ai]C = 0. (443) 

The subcase C = 0, h ence S = —602, ai — 0, yields the degeneracy ai = of 
the three fronts ( [43 2| ). In the other subcase, the system for C is linearizable 
into the third order system [S) ( [430| ) in which both d^ and dt change sign, 
the generic solution {S, C) of (_E2, -B3) is therefore 

bC = -3axLog!f3(-a;,-t), S = -602 + {bCf/Q-bC^, (444) 

and there only remains to integrate (|64|)-(|65|) for x or (|6C|)-(|6l|) for Lp. Since 
the one-form da; — Cdt possesses an integrating factor ^, the PDE ( |6l] ) for 
(/3 can be integrated by the method of characteristics, 

1'x+ k2^ ^-k2X~klt/b 

Note that the cyclic permutation of the roots Cj is broken when going from 
(5*, C) to tp. With the classical identity on Schwarzians 

{^■,x} = {1>-F}Fl + {F-x}, (446) 

the third order ODE (^) for ip becomes 

{<P;F}^Q, (447) 

which integrates as 

(p = (!>[F) = -^ -^, Aj arbitrary constants, AiA^ - A2A3 ^ 0. (448) 

A^iF + A4 

The value of x^^ 

F F 

X^^ = -Er^rr - T^- ^0 arbitrary constant, (449) 
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is again invariant under a cyclic permutation of the roots e, , and the solution 



u finally obtained is (431) 



The one- family truncation (142) with r = ■0 and the third order assump- 



tion (117)-(118) generates five determining equations (15C), their straight- 
forward resolution yields 



u = ddx Log?/; + si/3 -I- dU, 
'ipxxx + iUipxx - 3(a2 - [/^ - Ux)ll^x 
- {ai + 3a2U-U'^ -bUt/2- 
h^t - i^xx - 3f/i^, - 6L/V = 0, 



(450) 

WU^ + (7,,/2)V' = 0, (451) 
(452) 



in which U is constrained by two relations. But, since the coefficient / can be 
set to zero without loss of generality, the choice [7 = represents the general 
solution (just like i n |31[| where constant values were assumed ah initio for 
the coe ffici ents in (117)-(118)), and it represents again the collision of two 
fronts ( pl| ). 

The contrast of difficulty between the second order assumption (laborious) 
and the third order assumption (immediate) is the signature that the good 
scattering order of KPP is three, despite the irrelevance of such a notion for 
nonintegrable equations. 



The two- family truncation with a second order assumption (371)-(374) 
[P0pi,100| generates five determining equations. Despite the factorized form 
of £'5 , we have not yet found their general solution. The three particular so- 
lutions for which (S', C, A) are constant provide immediately the three pulses 
(435), for they belong to the class of polynomials in tanh and sech, generated 
by negative and positive powers of x according to the elementary identities 
(pTl). 



The elliptic solution ( |433| ) can also be written 

u = si/3 -1- dx Log(ns(V') - cs(?A)), V = !^3, 53 = 0, §2 arb., ai = 0, (453) 

a relation in which the argument of the logarithm is the ratio of two entire 
functions. Therefore it could be possible to find it by a suitable extension of 
a two-family truncation with a third order assumption. 

This solution was first found by the following two-step procedure pO| , 
which, unfortunately for this nice method, only works for a restricted class 
of PDEs (those with p = — 1, uq = cq, ui = ciC -I- C2, Cj = constant, see (|80|)). 
The first step is to define the truncation 



u — ddx Log((/3 — ifo) + U, U = constant. 



E{u) ^y^^Ej{ipxx/ipx,ipt/(px 
3=0 






i-3 



, Vj : Ej = 0, (454) 



whose general solution is U ~ si/3,(^— ipo — ^3 (indeed, comparing with 
(pOl), this assumption a priori implies bipt — Sipxx = 0). The second step is 
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not a truncation, but the change of function u t-^ f 

u = si/3+{dd^Log^3)fi^3), (455) 

which transforms ( jn\) into 

U" - 2[/3 + 2ai^.^^ = 0, Uiij) = /(^)/V. (456) 



This is an ODE iff ai = 0, in which case its solution is (433). Therefore, 
the assumption ( [455| ) has defined a reduction of the PDE to an ODE. This 
subject will be further examined in Section nG. 

The cubic complex Ginzburg-Landau equation 

The cubic complex Ginzburg-Landau equation (CGL3) 

E{u) = iut + pu.^j. + q\u'\^u — i^u = 0, pq ^ 0, {u,p,q) G C, j £ TZ, (457) 

with p, q, 7 constant, is a generic PDE describing the propagation of the sig- 
nal in an optical fiber as well as superfluidity, spatiotemporal intermittency, 
pattern formation, etc. 

One easily checks that |u| generically behaves like a simple pole. The 
dominant behaviour 

u^aox-'+*", ^~aoX~'~'", (458) 

in which ag is a complex constant, a a real constant, is solution of the non- 
linear algebraic system 

p{-l + ia){-2 + ia)+qa2 = 0, (459) 

p{~l-ia){~2~ia)+qa2^0, (460) 

with a2 — |aop. This defines two families for |up (four for |m|) ||^ 

q|„|2 o 

?-^dr- id,, A^ = dl + {8/9)dl (462) 

q 

To prevent these irrational expressions to mess up all subsequent computa- 



tions (Fuchs indices, no-log conditions, truncations), the system ( [459| )-(460) 
can equivalently be solved as a linear system on C [|3l| , p5| 

a2^--il-ia){2-ia), (463) 

q 

p = Kp{l - ia){2 - ia), q = Kq{l + ia){2 + ia), (464) 

in which K is an irrelevant arbitrary nonzero complex constant. 
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The indicial equation is the determinant 

P(j) ^ 

l + ia){j -2 + ia),p{j 



{2aoao)q 

■^ {2aoa^)q 



alq 



qao 
diag(p(j 



of the second order matrix 



1 ~ ia){j -2- ia)), (465) 



and with the resolution (463)- (464) it evaluates to 

detP(j) = (j + l)jif - 73 + Qo? + 12) 



0. 



For generic values of (p, q) , two of the four indices are irrational. 
Let us consider, for simplification, the solitary wave reduction 



u(:E,t) = C/(e)e<^^ + ^(^)) 



, t- 



ct. 



(466) 



(467) 



in which ([/, f) are functions of the reduced independent variable ^, and let us 
restrict to the pure CGL3 case Im(p/g) ^ 0. The general solution of the fourth 
order system of ODEs for ([/, ip) a priori depends on six arbitrary constants, 
the four constants of integration plus the two reduction parameters {c,ui). 
From these six constants, one must subtract 

1. the irrelevant origin ^o of ^ (Fuchs index —1), which represents the in- 
variance under a space translation, 

2. the irrelevant origin (fo of the phase (Fuchs index 0), which represents 
the invariance under a phase shift, 

3. and the number of irrational Fuchs indices, generically two. Indeed, these 
irrational indices represent the chaotic nature of CGL3 (see the expansion 
(48) in pl|) and they cannot contribute to any analytic solution. 

Therefore only two relevant arbitrary constants are present in what can be 



called the general analytic solution of the reduction (467). 

Presently, one only knows four particular solutions of the reduction ^ = 
X — ct with a zero codimension (no constraint on (p, g, 7)). These are 

1 . a pulse or solitary wave [M 



U — — iflQ 

2. a front or shock I 



fcsechfcxe* ["Log cosh fcx + i^ii]^ ^^2 



7 == 0, (468) 



""2 



tanh -^ ± 1 

2^ 



z [a Log cosh —^ 
e 2 



3. a source or propagating hole g 
{Ki + iK2)c 



flo 



k , k^ 
— tanh — f 
2 2^ 



i [a Log cosh — C + ^3C^ 
xe 2 



{Kik^ + K5C^)t] 



K^c^ - KicH] 



(469) 



Kek^ + Krc'^ =7,(470) 
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4. an unbounded solution [^ 

|u|2 = a2(tan2 -^ + K^), c = 0. (471) 

In the above expressions, all parameters {a2,a,k,c,Ki) are real and only 
depend on {p,q,j), except in ( |470| ) where the velocity c is arbitrary. 

These four particular solutions are four different degeneracies [pll of the 
yet unknown general analytic solution. 

In experiments or computer simulations, one has observed |5 y73| , |l4p6| 
other regular patterns which should correspond to other degeneracies of the 
general analytic solution. One of them p6| is a homoclinic hole solution, com- 



plementing the heteroclinic hole (470). Another o ne, o f the highest interest in 
fiber optics, is a propagating pulse, extrapolating (f46S| ) to c 7^ and reducing 
in the NLS limit {p, q real, 7 = 0) to a "bright soliton" of arbitrary velocity. 
Let us now address the question of retrieving these four solutions (and 
ideally of finding the unknown one ar at least other degeneracies) by some 
truncation. For a truncation to be successful, the truncated variables should 
be free of any multivaluedness in their dominant behaviour. This is not the 
case of the natural physical variables {u,u) or (Reu,Imu), which are always 



locally multivalued as seen from (458). A more detailed study [pl| uncovers 



the best representation for this purpose, namely a complex modulus Z and a 
real argument uniquely defined by 

u = Ze'^, u = Ze-'^, (472) 

and the above four exact solutions are written in this notation. For each 
family, if one excludes the contribution of the irrational Fuchs indices, the 
three fields {Z, Z, grad 0) are locally singlevalued and they behave like simple 
poles. The physical variables (|Mp,gradargu) also have this nice property of 
being locally singlevalued (they respectively behave like a double pole and a 
simple pole), but they are not as elementary as {Z, Z,gia,d0). 

The one- family truncation of the third order ODE satisfied by |up (after 
elimination of </?), with a constant coefficient second order assumption, ev- 
idently captures all four solutions, since jup is a degree- two polynomial in 
tanh K^. Such a truncation generates cumbersome computations and provides 
no additional solution. 

The one-family truncation of (Z, Z , grad6>) with the same constant coef- 
ficient second order assumption is defined as [pllJSq] 



(473) 



in which x and ^|: arc functions of ^ = x~ct, (a;, X, Y, K, k'^) are real constants. 
One has to solve the four complex (eight real) equations Ej = in the eight 



z ^ 


aoix 


-^+X + iY), 




z = 


aaix 


-^+X-iY), 




= 


-- ujt + 


aLogV" + K^, 




(Lo 


gV')' = 


'-x'\x'^i- 


{ky^)x' 


Ee- 


-i0 _ 


V"3 RyJ-3 
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real unknowns {a2,a,u>, X,Y, K,c,k'^), the two complex parameters {p,q), 
and the real parameter 7. If there exists a solution, the elementary building 
block functions evaluate to 

X = - tanh y , V = cosh y . (474) 

The good methodology |3|,^ is again to select, among the eleven complex 
variables considered as equivalent, four variables which make the system a 
linear one of Cramer type. The system {Eo,Ei,E2) is of Cramer type in 
{a2,K,(jj), and after its resolution the last equation E3 is independent of 
(p, g, 7, c) and factorizes into a product of linear factors 

E3 = [k^ - A{X + iYf] {aY - 2X) = 0. (475) 

Finally, this one-family truncation recovers all four solutions except the pulse 

(Eel). 

The two-family truncation of (Z, Z, giadO) with the same constant coef- 
ficient second order assumption retrieves the pulse solution (p6q), but finds 
nothing new. 

Similar truncations for two coupled CGL3 equations can be found in 
^5|M 



The nonintegrable Kundu-Eckhaus equation 

The PDE for the complex field U{x, t) ||,| 

iUt + aU^^ + {—\U\'' + 2he''^{\U\\)U^Q, ia,P,b,j)€n, (476) 
a 

with a/3 cos 7 ^ 0, is linearizablc when 6^ = /3^ into the Schrodinger equation 



iVt + aV,, = 0, [/ = ./ " ^!^ =. (477) 

V 2/3 cos 7 ^/j^^ 

This suggests considering the PDE for u — J \V\'^dx 

a 2 3 9 r,/5^ ^ (^sin7)^ 4 

TTK'^xxxx^x ' ^xx ^'^x'^xx'^xxx ) I ^ ^x^xx 

2 a 

+2{bcosj)uluor:xx + 7;-{uttul + u^xuf - 2utUxUxt) = 0. (478) 

2a 

When b^ 7^/3^, this PDE fails the test pSJ because, for each of the two families 
for u, one index is generically irrational. However, its one-family truncation 
with a second order assumption (i.e. the usual WTC truncation) is a very 
rich exercise [p2[ which yields quite unusual solutions, among them an elliptic 
one involving the ODE of class III of Chazy [n3| . 
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10 Singular manifold method versus reduction 
methods 

In order to find exact solutions of PDEs, there exist two main classes of 
methods. The first class, which has been detailed in these lectures, is based 
on the structure of the movable singularities and it can be called, in short, 
the singular manifold method. 



The second class, presented in another course of this school |121|, basically 
relies on group theory and consists in finding the reductions to a PDE in a 
lesser number of independent variables, and at the end to an ODE. These 
reductions are obtained either by looking for the infinitesimal symmetries 
of the PDE (space translation, etc) and by integrating them, or by a direct 
search not involving any group theory. The main methods in this second class 
are known as (see references in [ |l2l| ) the classical method (point symmetries), 
the nonclassical method (conditional point symmetries), the direct m,ethod 
(direct search). 

The question of the comparison of these four methods by their results is 
an active research subject Q,^ 52 , and its current state is given in ||l7| , |l2l| . 



Let us take as an example a second order nonintegrable PDE, this is 
enough to give an idea of the comparison. The KPP equatio n (|77| ) has been 
studied in detail with the singular manifold method. Section |9.3| . It has also 
been investigated with the three other methods, and the results are the fol- 
lowing. 

In the classical and nonclassical methods, let us denote 

t{x, t, u)ut + £,{x, t, u)ux — rj{x, t,u) = Q (479) 

the PDE for u{x,t) which, after computation of the symmetries (t, ^,77), 
defines the constraint on u susceptible to yield a reduction if the constraint 
can be integrated. 

Classical method . It yields only two reductions u{x,t) t— > U{z) P,[l0|, one 
noncharacteristic (i.e. conserving the differential order two) 

z = x-ct, u^U, -U" - bcU' + 2d-'^{U - ei){U - e2){U - 63) = 0, (480) 

one characteristic (i.e. lowering the differential order two) 

z^t, u = U, bU' + 2d-^{U - ei)([/ - e2)(C/ - eg) = 0. (481) 

This is in fact a unique reduction z — Xx + fit, but the splitting according to 
the characteristic nature is relevant for the Painleve property. None of these 



two ODEs has the Painleve property, unless c — {3ej — si)/{bd) in ( 480| ). 

Nonclassical method . It yields three sets of values for (r, ^,77) |96||, two 
with T 7^ and one with r = 0, ^ 7^ 0. 

The first one pU has codimension one 

ai = 0, ei = (62 + 63)72, but - 3d-2 ((LogV'),(u - 61))^ = 0, 1/; = <p3, (482) 
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its integration defines the noncharacteristic reduction to an elliptic equation 
ai = 0, z = 1^3, u = ei + dz^U{z), U" - 2U^ = 0, (483) 



and one finds the solution ( [433| ). 

The second one ^^ has codimension zero 

but + d^^{u - si/3)uj: + 3d"2(u - ei)(u - e2)(w - 63) = 0, (484) 

and, remarkably, this first order PDE, which fails the t est, identifies to the 
no-log condition (82) with 3m— si — bdC. Its integration ( 436 ) cannot define a 
reduction unless some choice of the arbitrary function is made. Nevertheless, 
the common solution to the PDEs (^) and ( |484[ ) is ( |43l[ ). 
The third one M is 



Ux - ri = 0, 
in which 77 satisfies the second order PDE 

r]xx + '^V'llxu + Tfvuu + 2d^^{u - ei)(u - e2){u - e3)riu 
-2d-2(3y2 _ 2s^y + 52^3 _ 3(fa2)r] -\- brjt = 0. 



(485) 



Integrating ( 486 ) is equivalent to integrating the original PDE (|77|), since the 
transformation ( |485| ) simply exchanges them, so one is stranded. The only 
way out is to put some additional constraints on 77. The consistent way to do 
that (p7[ page 634) is to eliminate u^ and its derivatives (in this case u^x 
only) between ( [485| ) and (^), which results in a nonlinear first order ODE 
for the function t h^ u{x, t) (i.e. with a; as a parameter) 



but - iVx + VVu) + 2d ^{u - ei)(u - e2)(u - 63) = 0. 



(487) 



Requiring the invariance of this ODE under the infinitesimal transformation 
(r = 0,^ = 1,77) of the classical method creates constraints on 77, an exercise 
which is left to the reader. 

Direct method . The search for a reduction u{x,t) 1-^ U{z) in the class 

u[x, i) = A(a;, t)U{z{x, t)) + ^(a;, i), (488) 

apart from the characteristic reduction z = t,u = U , yields the noncharac- 
teristic reduction M 



U"-2U^+g,{z)U' + g2{z)U + g3(z)^0, 
u — dzxU + si/3, Zx 7^ 0, 

provided 2,171,(72,53 satisfy the system 

- Qa2Zx + zlg2 = 0, 

+ zlgi = 0, 

= 0. 




(489) 
(490) 



(491) 
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The constraint Zj; 7^ splits the discussion into oi = and ai 7^ 0. The case 
oi arbitrary defines the reduction p& 

z = x-ct, u = dU + si/3, U" + cU' - 2{U^ - 3a2U - ai) = 0, (492) 



identical to (48C). In the case ai = 0, hence 53 — 0, the system is solved for 
(zj^xxyZt), and the condition {zxxx)t = {zt)xxx reads 



18 



18- 



+ 2g[{z) - gi{z)d, + 39^ \q{z)^ 0, (493) 



Q(z) = 9.92-2.g2-35i. 
For Q{z) ^ 0, the condition integrates as 

Z = G(x/i(t)+/2(t)), 



(494) 



(495) 



in which G is an arbitrary function, and (/i,/2) arc further constrained by 
/{ = 0, /^ = 0. The resuh is ai = 0, z = G{x), and the ODE ( |489|) transforms 
to an elhptic equation under 



C/(^) - ^^^§7^, /"-2f + 6a2/ = 0, 



(496) 



in which G ^ denotes the inverse function of G. This solution is not distinct 



from the stationary elliptic reduction (434). 

For Q{z) — 0, if one defines the function G by 

5i(z) = -3(Log(G'(z)))', 

the system (|49l[) is equivalent to 



fff2 = (2/9)g? + (l/3)g'i, 

{dl--ia2d^)G{z{x,t))^(), 

{bdt-idl)G{z{x,t))^0, 
[ f"{Z) - 2f{Zf = 0, U{z) = G'{z)f{Z), Z = G{z), 



(497) 



(498) 



and this proves that the particular solution gi — g2 — ^ considered in [ p6[ is 
the general solution, equivalent to the reduction z = if's in (433). 



11 Truncation of the unknown, not of the equation 

When applied for instance to the second Painleve equation (P2) 

(P2) E{u) = u" -2u^ -xu-a^O, (499) 

the one-family singular manifold method in the case of a second order s cat- 
tering problem, i.e. the one originally performed by WTC, see Sections 7.1 
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and 7.3, presents the following drawback. The (P2) ODE has two families 
u ~ ex~^, e^ = 1- With the definition of gradx 

X' = 1 + (5/2)x^ (500) 

the one-family truncated expansion for u is found to be 

u = ex-\ E{u) = e{S - x)x-^ + (-a - eS' /2)x\ (501) 

E2 = e{S - x) = 0, (502) 

E3, = -a- eS"/2 = 0, (503) 

and its general solution is 

S = x, 2a + £ = 0, u = e(Log Ai(x))', Ai" +(.x/2) Ai == 0. (504) 

One therefore finds only a one-parameter particular solution in terms of the 
Airy function, at the price of one constraint on the parameter a. This is 
unsatisfactory because the method fails to find the highest information on 
(P2) (highest in the context of these lectures), namely its Schlesinger trans- 
formation. Such a transformation is by definition a birational transformation 
between two different copies of (P2), denoted u{x,a) and U{X,A), and it 
reads [fl 

X = X, u + U ^ — ^- = — TTs , a = A + 1. (505) 

2{U' + U^)+X 2(u'-u2)-a;' ^ ' 

A method to remedy this drawback is the following |lq|. We rephrase it 
in the homographically invariant formalism, which simplifies the exposition. 
Firstly, rather than splitting E{u), defined in (p01[), into one equation per 
power of Xi one retains the single information E(u) — 0, and one eliminates 



u and X between the three equations ( |50C| ) and (501) to obtain the second 
order ODE for S{x) 

2{S - x)S" - S"2 + 25" + 2S^ - 4xS^ + 2x^S + 4a(a + e) = 0. (506) 

This ODE for S(x), which is birationally equivalent to (P2) under the trans- 
formation 

u = ex-\ S = -2(x-^)' - 2x-^ x'' = ^^, (507) 

2(5 — x) 

bears the number 34 in the classification of Painleve and Gambler |£l| . 

Secondly, despite the fact that one already knows the general solution S{x) 
in terms of the (P2) function u(x^a), one takes advantage of the two- family 
structure of (P2) (the sign e is ±1) to perform an involution by representing 
S{x) with another (P2) function U{X,A) as 

S = -2V' ~ 2t/^ V = eaV", e\ = 1, V" -2U^ -XU - A = Q, X = a;.(508) 
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The elimination of S between (507) and (508) provides a relation between 
{ea,e2A) only 



{A + e2/2f = ia + e/2f 



(509) 



The solution A = —e2{sa + 1) is the Schlesinger transformation. 

An equivalent presentation can be found in Ref. p5| . In the latter, one 
first computes the two coefficients uo,ui of the Laurent expansion 



u = uox ^ +ui, 



(510) 



then the Schlesinger transformation is readily obtained by (more precisely, 
the computation of |^ reduces to) the elimination of the three variables 
u, Z, U" between the four equations (mq, mi, u, U, Z, s are functions of X = a;) 



u = uqZ ^ + U, 



Z' = \ 



U-ui 



-Z\ 



(511) 

(512) 

(513) 
(514) 



uo 2 

(Pn)(u,a;,a,/3,7,5) = 0, 
(Pn)(C/,X, A, B,C,i:') = 0. 



Equation ( pllj) is an assumption for a Darboux transformation, and ( pi 21 ) 
defines a Riccati equation for the expansion variable Z which depends on a 
free function s. The elimination is differential for u and Z, algebraic for [/", 
and it results in 



F(C/', f/; s, s', s", a, /3, 7, (5, A, B, C, D) = 0. 



(515) 



The algebraic independence of {U',U), consequence of the irreducibility of 
(Pn), requires the identical vanishing of _F as a polynomial of the two variables 
{W, U), and this provides two solutions: the identity {u = U, Z^^ = 0) and, 
at least for (P2) and (P4), the Schlesinger transformation. The result for (P2) 
is 



(P2) eZ-^ =u-U = 



e{A- a) 



2U' + e(2C/2 + x) ' 
and the inverse transformation 

e{A-a) 



(P2) u-U 



2u' + e(2u2 + x) 



^ + £ = 0, s = 0, (516) 



(517) 



follows from the elimination of U' between (516) and 

(L/-w)' + e(C/2-M2) = 0, 



(518) 



itself obtained by the elimination of Z between (511) and (512). This Schlesin- 
ger t ransformation is identical, thanks to the parity invariance of (P2), to 

( Hosl) . 
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The result for (P4) is 
(P4) u" ~ u'V(2u) - (3/2)^3 - ixu^ - 2x^u + 2au - P/u = 0, (519) 



3U' + e(3C/2 + 6xU - 2^ - 4a) + 6 ' 

4£(a - A)u 



~ 3w' + e(3u2 + 6xu - 2a - AA) + 6' ^ ^ 

{U - u)' + e(t/2 -u^ + 2x{U -u) + 2{a - A)/3) = 0, (521) 

9/3 + 2(a + 2A - 3e)2 = 0, 9B + 2{A + 2a- Sef = 0, (522) 

s = 4(A-a)/3. (523) 

12 Conclusion, open problems 

The singular manifold method, which is based on the singularity structure, 
is quite powerful to provide exact solutions or other analytic results. There 
still exist many challenging problems, in particular in nonlinear optics and 
spatiotemporal intermittency |5y56|] , in which the equations, although nonin- 
tegrable, possess some regular "patterns" which could well be described by 
exact particular solutions. The difficulty to find them |Q comes from the 
good guess which must be made for the functions ip, which do not necessarily 
satisfy a linear system any more. Methods from group theory usually pro- 
vide complementary results, although they also fail in the two just quoted 
examples. 
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